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HYPOXANTHINE, A GROWTH SUBSTANCE FOR PHYCOMYCES 
By WILLIAM J. RoeBINS AND FREDERICK KAVANAGH 
NEw YorRK BOTANICAL GARDEN AND DEPARTMENT OF BoTaNy, COLUMBIA UNIVERSITY 
Communicated January 23, 1942 


In previous papers from this laboratory an unidentified growth sub- 
stance for Phycomyces widely distributed in products of natural origin was 
described as factor Z;.!_ Guanine was found to have the same effect as 
factor Z;, but was not considered to be identical with it because the activity 
of the former was destroyed by treatment with nitrous acid while that of 
the latter was unaffected. Xanthine, as might be anticipated from the 
effect of nitrous acid on the activity of guanine, was found to be inactive. 
Choline, adenine, thymine, uracil and cytosine also were ineffective.? 
Further experiments have been carried out in an attempt to identify fac- 
tor 2}. 

The procedure followed was that described in the previous paper.” The 
growth of Phycomyces during a period of 72 hours at 26°C. was observed in 
25 ml. of a basal mineral-dextrose solution containing asparagine and thi- 
amine to which the various compounds or preparations tested for factor Z; 
activity were added alone and in combination with a Dr fraction* prepared 
from white potato tubers. We have used also the 48-hour development in 
10 ml. of solution, a modification which economizes time and material. 

The following were tested for factor Z, activity: monomethyl guanosine, 
monoacetone guanosine, guanosine picrate, guanylic acid, sodium guanyl- 
ate, isoguanine sulfate, sodium inosine, hypoxanthine, vicine, divicine sul- 
fate, barium inosinate and 2-amino uric acid (all obtained through the 
courtesy of the Rockefeller Institute for Medical Research), hypoxanthine 
(prepared in our laboratory by treating adenine with nitrous acid), a con- 
centrate of folic acid supplied by R. J. Williams and one obtained from W. 
H. Peterson, uric acid, amygdalin, allantoin, theobromine, theophyllin, 
theocin,‘ adenine, guanine and xanthine. 

Guanine and both samples of hypoxanthine were active and their ac- 
tivity on an equimolecular basis was approximately the same (table 1). 
None of the other compounds or the two concentrates of folic acid was 
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effective, at least in amounts of the same order of magnitude as the effective 
quantities of guanine or hypoxanthine. Hypoxanthine® was isolated in 
crystalline form as the silver salt (Fig. 1) from a Ca fraction from white 
potato tubers which had been treated with nitrous acid and electrodialyzed 
and which showed factor Z; activity. From 60 ml. of the solution 15.5 mg. 
of the silver salt of hypoxanthine were obtained. The substance isolated 
from this salt was active. On this basis it was estimated that the factor Z; 
activity of the extract from which the hypoxanthine was isolated could be 
accounted for by its hypoxanthine content. It would appear, therefore, 
that the active substance present in the potato extract after treatment 
with nitrous acid was hypoxanthine. 


TABLE 1 

THE EFFECT OF GUANINE AND HYPOXANTHINE IN THE PRESENCE OF THIAMINE AND A 

Dp FRACTION FROM Potato TUBERS. Notice THAT BoTH ComMPOoUNDS ARE ABOUT 
EQUALLY EFFECTIVE IN INCREASING THE BENEFICIAL ACTION OF THE Dp FRACTION 


ADDITIONS TO 10 ML. DRY WT. ADDITIONS TO 10 ML. DRY WT. 


OF BASAL SOLUTION 2 CULTURES OF BASAL SOLUTION 2 CULTURES 
PLUS TRTAMINE MG, PLUS THIAMINE Mc, 
10 ug. guanine HCl 6.6 ug. hypoxanthine 
plus 63 plus 64 
Dp fraction Dp fraction 
1 yg. guanine HCl 0.66 ug. hypoxanthine 
plus 61 plus 67 
Dp fraction Dp fraction 
0.1 wg. guanine HCl 0.066 ug. hypoxanthine 
plus 43 plus 42 
Dp fraction Dp» fraction 
None 0.5 Dz fraction 27 


The active material present in the original potato extract also appeared 
to be hypoxanthine. Examination of an active Ca fraction showed little 
or no guanine or xanthine, some adenine and considerable hypoxanthine.® 
The quantity of hypoxanthine found in this extract was sufficient to account 
for the activity of the fraction. The hypoxanthine crystallized from the 
potato extract was active; the adenine was inactive. Furthermore, the 
Ca fraction was as active (or more active) after treatment with nitrous acid 
as before. This would be expected if the active material in the original 
potato extract were hypoxanthine. We are inclined to believe, therefore, 
that factor Z; in potato extract is hypoxanthine though in other materials 
it may be guanine or a mixture of the two. 

Hydrolysis of guanosine, guanylic acid, sodium guanylate, guanosine 
picrate, monomethyl guanosine, sodium inosine, barium inosinate and 
monoacetone guanosine yielded active preparations probably because of 
the guanine or hypoxanthine formed in the hydrolysis. Both samples of 
folic acid also gave active material on hydrolysis, that from Peterson being 
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the more active. This indicates that folic acid probably contains guanine 
or hypoxanthine. 

The small quantity of hypoxanthine which is effective (activity has been 
observed for 0.3 my mole) is evidence for its functioning as a growth sub- 
stance. Furthermore, the inactivity of xanthine, adenine and other com- 
pounds closely related to guanine or hypoxanthine demonstrates a consider- 
able degree of specificity for the action of these two compounds and sug- 
gests that they may function in an enzyme system as thiamine, riboflavin, 
adenine and nicotinic acid are known to function. The inactivity of 
xanthine shows that the replacement of the NH, radical or hydrogen in the 
second position on the purine ring with oxygen renders guanine or hypo- 
xanthine inactive (table 2). The inactivity of adenine suggests that oxy- 





FIGURE 1 


(1) Crystals of the silver salt of hypoxanthine prepared from pure hypoxanthine. 
(2) Crystals obtained with ammoniacal silver nitrate from potato extract which had 
been treated with nitrous acid and which showed factor Z; activity. Magnified 160 x. 


gen in the sixth position is important because the substitution of NH2 for 
oxygen in hypoxanthine renders that compound inactive. The inactivity 
of 2-amino uric acid suggests that the hydrogen in the eighth position is 
important and may not be replaced by oxygen. The inactivity of guanylic 
acid, guanosine and similar compounds indicates that position 7 or 97 must 
be open, perhaps for combination with some other compound than the 
sugar which is present in guanosine. The inactivity of isoguanine empha- 
sizes the importance of oxygen in the sixth position and hydrogen or NH: 
in the second position. It would appear that the active substance must 
have the eighth position and the seventh or ninth position open for sub- 
stitution, oxygen in the sixth position and hydrogen or NH, in the second 
position. : 
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We can only speculate as to whether hypoxanthine (or guanine) is active 
per se, or whether it, together with an unidentified substance, makes up a 
larger and active molecule, as, for example, thiazole and pyrimidine com- 
bine to form thiamine. Our observations have shown that factor Z; 
(hypoxanthine or guanine) is much more effective in the presence of a sec- 
ond unidentified substance, factor Z2, than it is when used alone. Further 
information on the relation between factor Z, and factor Z, will have to 
wait the isolation of the latter substance. 


TABLE 2 


POSITION ON THE PURINE RING OF VARIOUS RADICALS FOR SOME OF THE COMPOUNDS 
TESTED FOR Factor Z; ACTIVITY 
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SUBSTANCE POSITION ON PURINE RING ACTIVITY 
1 2 3 6 7 8 9 
Guanine H NH: ee O H H Active 
Hypoxanthine H H O H H Active 
Isoguanine = O H NH: H H Inactive 
Xanthine H O H O H H Inactive 
Adenine H i NH; H H Inactive 
2-amino- 
uric acid H NH: 6A O H O H Inactive 
Theobromine H O CH; O CH; H Inactive 
Theophyllin CH; oO CH; O H H Inactive 
Guanosine H NH: ni O sugar H Inactive 


1 Robbins, W. J., Am. Jour. Bot., 26, 772-778 (1939); Robbins, W. J., Bot. Gaz., 101, 
428-449 (1939); Robbins, W. J., Am. Jour. Bot., 27, 559-564 (1940); Robbins, W. J., 
and Hamner, K. C., Bot. Gaz., 101, 912-927 (1940); Robbins, W. J., Zbid., 102, 520-535 
(1940). 

2 Robbins, W. J., and Kavanagh, F., Proc. Nat. Acad. Sci., 28, 4-8 (1942). 

3 The Ca fraction was the material adsorbed from an extract of potato tubers on 
charcoal and eluted with ammoniacal acetone. The Dz fraction was the filtrate from 
the charcoal treated extract. 

4 Theocin is a trade name for theophyllin. Our preparation of theocin was synthetic. 

5 The crystals of the silver salt of hypoxanthine were prepared as follows: A Ca frac- 
tion from potato tubers was treated with nitrous acid and electrodialyzed. The cathode 
fraction was evaporated to small volume and the purines precipitated with ammoniacal 
silver nitrate. The precipitate was dissolved in dilute hot nitric acid. The solution 
was cooled, the precipitate removed and redissolved in hot dilute nitric acid. This 
procedure was repeated five times. , 

6 The precipitate formed from 500 ml. of a Ca fraction by ammoniacal AgNO; was 
freed from silver by treatment with HCl. Picric acid was added and the precipitate 
filtered off. Ammoniacal AgNO; was added to the filtrate and the precipitate was dis- 
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solved in hot dilute HNO;. The precipitation with AgNO; and solution in HNO; was 
repeated five times. The pale yellow solution was treated with a small amount of Norit 
A and filtered. From the filtrate 15 mg. of crystals were obtained identical in appear- 
ance with the silver salt of hypoxanthine. 

The precipitate of picrates was dissolved in dilute NH,OH and treated with AgNO. 
The precipitate was filtered off, washed, decomposed with HCl to remove the silver and 
filtered. The filtrate was made slightly alkaline with NH,OH and allowed to stand for 
24 hours. No guanine precipitated. It was treated with nitrous acid to convert any 
guanine and adenine to xanthine and hypoxanthine, respectively. Silver nitrate was 
added and the precipitate examined. It appeared to consist of the silver salt of hypo- 
xanthine. In the examination of the picrate precipitate volumes were kept to less 
than 20 ml. 

™ There appears to be some uncertainty as to whether the sugar in guanosine is at- 
tached to the purine ring in the seventh or the ninth position. It is usually represented 
in the seventh position. 


SIZE GENES OF MICE 
By W. E. CastTLe 
DEPARTMENT OF GENETICS, UNIVERSITY OF CALIFORNIA 
Communicated January 21, 1942 


In previous papers! it has been shown that several mutant genes of the 
house mouse, in addition to other effects which they exert, act also as size 
genes, either increasing or decreasing general body size. 

The present paper contains a report on a study of the effects of two coat 
color genes, pink eyez (f2) and lethal yellow (A”) on body size, separately 
and in association with each other, but always on a background of homo- 
zygous brown (bb). Ina previous paper it was shown that the influence of 
p2 on body size is increased by association with brown. The same is proba- 
bly true of A” also, although qualitatively its action is the opposite of that 
of po, since it increases rather than decreases body size. 

In order to learn more about the interaction of the three mutant genes, 
pe, A” and b, crosses were arranged which would produce as litter mates 
mice of the four genetic classes (1) pink-eyed brown (aa bb pope), (2) pink- 
eyed yellow (A” a bb pops), (3) brown (aa bb Pe pe) and (4) yellow (A” a bb 
Pe p2). 

All four groups would be homozygous for 6. Groups (1) and (2) would 
be homozygous for ~2 also, but groups (3) and (4) heterozygous for pp. 
Groups (1) and (3) would be free from the A” gene, but (2) and (4) would 
contain it. All possible combinations of 2 and A” would thus be obtained 
associated uniformly with bb. The symbol for brown will accordingly be 
omitted from the genetic formulae hereafter. 
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In this study, as in earlier ones, the young mice were separated as to sex 
at weaning, and kept until six months old, when they were weighed, chloro- 
formed and measured as to body length and tail length. 

In table 1 is shown for each sex separately the average size of the mice of 
each group as indicated by the three criteria, weight, body length and tail 
length. 

Group (3) consisting of aa P2: p2 animals shows body size influenced pre- 
sumably by the brown gene alone of those under consideration, since the A” 
gene is absent and the p2 gene present only as an unseen recessive, probably 
with little if any influence on body size. This group then should give us a 
standard for comparison of the effects of the genes p. and A” in various 
combinations. Average body weight of males in group (3) is 31.98 g., 
average body length 96.26 mm., average tail length 83.00 mm. The data 
for group (3) females are less satisfactory because of insufficient numbers. 
For the present we may pass them by. 

Group (1) consisting of aa 2 p2 animals shows body size influenced by 
both the brown gene and the pink eye gene in association. Average weight | 
of males in this group is 29.41 g., average body length 94.76 mm., average 
tail length 81.40 mm. These are decreases (from the group (3) standard 
averages) of 8 per cent in weight, 1.5 per cent in body length and 2.0 per 
cent in tail length, which must be ascribed to the influence of 2 now homo- 
zygous. Decrease is shown by females also of group (1) as regards body 
length (0.2 per cent) and tail length (2.9 per cent) but not as regards weight. 
Less importance attaches to these differences because of the insufficient 
number of females in group (8). 

Group (2) consisting of A” a pe» p2 individuals differs from group (1) only 
in the substitution of an A” gene fora. A marked increase in body size 
results amounting in males to 20.9 per cent in weight, 2.2 per cent in body 
length and 1.3 per cent in tail length. In females the corresponding in- 
creases are 23.2 per cent in weight, 2.3 per cent in body length and 1.4 per 
cent in tail length. The A” gene, as is well known, increases weight chiefly 
by leading to accumulation of fat especially in females, hence the large 
percentage of increase in weight. 

Group (4), consisting of A” a P» ps individuals, gives us by comparison 
with group (3) another opportunity to measure the influence of A” in sub- 
stitution fora. In both these groups the influence of 2 is negligible, since 
it is present only as an unseen recessive. The differential factor is A’. 

Group (4) males show an increase over group (3) males of 16.9 per cent in 
weight, 1.9 per cent in body length and 2.8 per cent in tail length. For 
females the corresponding increases are 38.1 per cent in weight, 4.4 per cent 
in body length and 1.9 per cent in tail length. But these figures are less 
reliable because of the small number of females in group (3) as already 
noted. 
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The relative body size of mice of each of the four groups, taking group (3) 
animals as a standard of comparison, is shown in table 1. Also by com- 
bining the data on groups (1) and (2) we obtain average values for all ani- 
mals homozygous for ~2, which may be compared with the averages for all 
animals heterozygous for p2, obtained by combining the data for groups (3) 
and (4). These combination values are shown in the lower half of table 1. 

By combining in a similar way the data for groups (1) and (3), we obtain 
values for all animals free from the A” gene (aa animals) which may be com- 
pared with the values obtained by combining groups (2) and (4), which in- 
clude all animals possessing the A” gene. See the lowest section of table 1. 

We thus obtain a comparison, based on all the available data, of the 
effectiveness of p2 in decreasing body size and of A” in increasing it, from 
which it appears that the influence of A” is greater on weight and body 
length but not on tail length, which, however, is least valuable of the three 
criteria of body size. 

The net percentage changes effected by the two genes are as follows: 


IN WEIGHT IN BODY LENGTH IN TAIL LENGTH 
by pe pz — 6.0% —1.4% —2.9% 
by A” a +24.7% +2.7% +1.8% 
Difference +18.7% +1.3% —-1.1% 


These figures obtained by summarizing all the data furnished by the four 
groups agree well with those based on group (2) alone, in which both pz pe 
and A” a were simultaneously operative, which are for weight a net per- 
centage change of +17.7, for body length +1.3 and for tail length —1.1. 

The results obtained in this experiment are similar both qualitatively 
and quantitatively to those obtained from earlier experiments. They show 
that the gene 2. when homozygous decreases body size as judged by the 
three criteria, weight, body length and tail length. They show also that 
the gene A” when heterozygous increases body size, notably in weight 
because of “‘adiposity,”’ but also though in lesser degree in body length and 
tail length. When both influences are simultaneously in operation, the 
influence of A” surpasses that of p: f2, much as regards weight, slightly as 
regards body length, but little or none as regards tail length. 

Summary.—The mutant gene, pink-eyez, decreases general body size in 
mice less than lethal yellow increases it, as regards both weight and body 
length, but not as regards tail length. 


1 See Castle, W. E., Genetics, 26, 177-191 (1941). 
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THE OXIDATION OF 3,5-DIIODOT YROSINE TO THYROXINE* 
By TREAT B. JOHNSON AND LYNDON B. TEWKESBURY, JR.T 
DEPARTMENT OF CHEMISTRY, YALE UNIVERSITY 


Communicated February 2, 1942 


A physiologically active protein product from which thyroxine (IV) can 
be isolated was prepared from milk casein by W. Ludwig and P. von 
Mutzenbecher! in 1939, by treatment of this protein with iodine under 
carefully controlled conditions. The experimental technique of these 
workers was repeated by C. R. Harington and Rosalind V. P. Rivers? in 
England, and the original findings of the German investigators were con- 
firmed in every respect. There appears to be no doubt as to the correctness 
of the conclusions reached by these independent workers and it, therefore, 
becomes a matter of immediate chemical interest to consider any possible 
mechanism of reaction whereby such an important biochemical trans- 
formation can theoretically be brought about. 

That iodine under proper experimental conditions is capable of affecting 
an oxidative coupling of two molecules of 3,5-diiodotyrosine to form a 
thyroxine molecule was established later by P. von Mutzenbecher,* who 
accomplished this change successfully by oxidation of 3,5-diiodotyrosine 
with hypoiodous acid formed during prolonged incubation of this a-amino- 
acid in alkaline solution at 37°. These experiments were repeated and the 
results reported by P. von Mutzenbecher?® were confirmed by P. Block, Jr.,* 
in 1940. The original hypothesis advanced by C. R. Harington? that the 
biogenic synthesis of thyroxine probably results from the iodination of 
tyrosine functioning in ‘the protein molecule through the intermediate 
stage of 3,5-diiodotyrosine must, therefore, receive careful consideration in 
the light of these most interesting results. 


Biogenic Synthesis: 
Tyrosine ——— 3,5-Diiodotyrosine ———> Thyroxine. 


The authors have repeated the incubation experiments of P. von Mutzen- 
becher® using synthetical 3,5-diiodotyrosine of the highest purity, and have 
confirmed his conclusion that thyroxine is a product of reaction under 
specific experimental conditions. If the postulation be correct that thy- 
roxine is produced here as a result of oxidation of the 3,5-diiodotyrosine by 
hypoiodous acid (HIO) formed from this a-amino acid in the alkaline solu- 
tion, then the yield of thyroxine should theoretically vary according to the 
quantity of active hypoiodate present in the solution. The authors, there- 
fore, modified P. von Mutzenbecher’s technique* by adding hypoiodate to 
the a-aminoacid solution in molecular proportion with the result that the 
yield of thyroxine was slightly increased as postulated. 
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After applying several oxidation experiments with successful results, and 
after giving more attention to the nature of the secondary products of the 
reaction, the authors have been able to formulate a reaction mechanism 
which has not hitherto been presented to explain the formation of thyroxine 
(IV) from 3,5-diiodotyrosine (I). Our scheme of transformation is based 
on the results obtained in the fundamental studies of R. Pummerer® and co- 
workers dealing with the manner of oxidation of o- and p-substituted 
phenols in alkaline solution. We submit, therefore, the following graphical 
representation of this reaction (table 1) leading to the formation of thyrox- 
ine (IV) from 3,5-diiodotyrosine (I). 


TABLE 1 
I 
2 HOC _CHrCH(NH,)COOH —2H 
I 
I 


o 


CH.-CH(NH;)COOH CH» CH(NH:)COOH 


II 
O 
| 
1 \I 
I 
CH, O0€ CH CH(NH,)COOH. etiaetay 
| 
CH(NH,)! 
COOH 
Il 
I I 
HO< >0€ CH: CH(NH,)COOH + 
i I 
IV 


CH, = C(NH,)COOH ———> CH,C(:NH)COOH #20 


CH;CO-COOH + NH; 
V 


A comparison of the graphical representation of R. Pummerer’s®* de- 
hydrogenation of p-cresol (VI) by oxidation with KsFe(CN), at 0° (table 2) 
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furnishes significant data which finds application in support of the hy- 
pothesis proposed in table 1. 





TABLE 2 
O 
OH T Setiks 
5 —2H - sesame 
——> 
CH; CHs 
VI 
O 


Rearrangement | 
———— HH 


cH; o> CH: CH; 0€ cH: 


VIII IX 


The chief product of this reaction is the tetrahydro-benzfurane structure 
expressed by formula (IX), and results apparently by rearrangement of 
the primary product of oxidation of p-cresol (VI) namely, the p-quinolether 
CysHyO-2 (VIII). In fact, according to R. Pummerer it is characteristic of 
o- and p-alkylated phenols to undergo dehydrogenation easily in alkaline 
media to form free radicals which dimerize to more or less stable aromatic 
quinolethers. Such free radicals are expressed by formulae (II) and (VII), 
and the dimerized products by formulae (III) and (VIII) in tables 1 and 
2, respectively. That o- and p-alkylated phenols may react in such a 
tautomeric manner to give derivatives of dihydrobenzene was also shown 
earlier by K. von Auwers® in 1902, who synthesized the quinone structure 
(X), for example, by interaction of p-cresol with chloroform in alkaline 
solution; and also by T. Zincke’ who demonstrated the production of the 
corresponding 1-methyl-1-trichlormethyl-4-ketodihydrobenzene (XI) by 
application of a Friedel and Craft’s reaction with p-cresol and carbontetra- 
chloride. 


0 O 
| | 
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“\ 
CH, CHCl CH, CCh 
x xI 
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If the hypothetical, quinolether intermediate (III), represented as 
formed from 3,5-diiodotyrosine (I) (table 1), is compared structurally with 
the corresponding oxidation product of p-cresol (VIII) (table 2), it may be 
seen that the internuclear condensation or rearrangement, which takes 
place in the formation of the tetrahydro-benzofurane (IX) from the p- 
quinolether (VIII), is precluded in the case of the corresponding hypo- 
thetical 3,5-diiodotyrosine derivative (III). The substitution of iodine 
atoms in both nuclear positions ortho to the ether-bound oxygen atom in 
(III) prevents any possibility of cyclization of this compound. Hence, 
stabilization of the system expressed by structure (III) must follow one of 
two courses namely, (1) molecular dissociation with loss of one alanine 
side chain and formation of thyroxine (IV) and iminopyruvic acid or, (2) 
hydrolysis of (III) with production of serine, HOCH2,CH(NH2)COOH and 
thyroxine (IV). The results obtained by experimentation to date lead the 
authors to the conclusion that the formation of thyroxine (IV) from 3,5- 
diiodotyrosine (I) as described above is an oxidation process, and that the 
correctness of the reaction-mechanism postulated is supported by the facts 
that pyruvic acid® (V) and ammonia have been identified by the authors as 
secondary products of the reaction. Thus far, the authors have been un- 
able to detect the presence of the a-aminoacid-serine as a product of the 
oxidation reactions. A thorough search for other products of reaction is 
now in progress. 

Summary.—It is a well-known fact that thyroxine and also 3,5-diiodo- 
tyrosine slowly undergo decomposition especially when exposed to sunlight, 
in weakly alkaline solution with formation of hypoiodous acid (HIO).® 
This halogen acid interacts slowly with 3,5-diiodotyrosine to form thyrox- 
ine,* 4 but thus far no mechanism of reaction has been approved to explain 
this biochemical change. The authors have repeated the work of previous 
investigators of this important change, and have found that not only 
thyroxine but also pyruvic acid and ammonia are formed as secondary 
products of reaction. The discovery of pyruvic acid (V) as a product of this 
oxidation reaction has made possible the postulation of a simple and plau- 
sible mechanism to explain the formation of the thyroxine. The study of this 
problem is being continued and a discussion of the experimental technique 
and analytical results will be published elsewhere. 


* Work supported in part by a grant from the George E. Sheffield Research Fund of 
the Sheffield Scientific School of Yale University. 

{ Present address: 222 Summer Street, Boston, Massachusetts. Ph.D. from Yale 
University in June, 1941. 

1 Ludwig, W., and von Mutzenbecher, P., Zeit. physiol. Chem., 258, 195-211 (1939). 

2 Harington, C. R., and Rivers, R. V. P., Nature, 144, 205 (1939). 

3 yon Mutzenbecher, P., Zeit. physiol. Chem., 261, 253-256 (1939). 

4 Block, P., Jr., Jour. Biol. Chem., 135, 51-52 (1940). 
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7 Zincke, T., and Suhl, R. S., Zbid., 39, 4148-4153 (1906). 

8 The pyruvic acid is easily identified by the successful application of A. Baeyer’s 
reaction involving the action of o-nitrobenzaldehyde on pyruvic acid which leads to the 
production of indigo (Baeyer, A., Ibid., 15, 2856 (1882); for other applications of this 
indigo color test see Johnson, T. B., and Baudisch, O., Jour. Amer. Chem. Soc., 43, 2670- 
2675 (1921); Ber., 55, 18-21 (1922). 
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COMPLEMENT FIXATION WITH SIMPLE SUBSTANCES 
CONTAINING TWO OR MORE HAPTENIC GROUPS 


By DAvip PRESSMAN, DAN H. CAMPBELL AND LINUS PAULING 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY, CALIFORNIA INSTITUTE OF TECH- 
NOLOGY 


Communicated February 13, 1942 


It has been shown that simple compounds containing two or more 
haptenic groups undergo: typical precipitin reactions with the antibody to 
a protein antigen containing homologous haptenic groups.' We have in- 
vestigated the question of whether or not complement fixation takes place 
in the formation of precipitates with these simple compounds and have 
obtained the results reported below. 

The compounds studied are listed in table 1. Phenylarsonic acid was 
included, although it does not form a precipitate with antiserum, because 
it is known to combine with antibody and thus to inhibit the formation of 
precipitate by the antibody and an antigen containing two or more benzene 
arsonic acid groups. 

Experimental.—The substances I to IV and the rabbit antiserum homolo- 
gous to azo sheep serum made with arsanilic acid were prepared by methods 
described elsewhere.'° The complement of the antiserum was inactivated 
by heating at 57° for 30 minutes. 

The amboceptor was prepared (by Mr. Carol Ikeda) by giving rabbits 
ten injections of 1 or 2 ml. of a 10% sheep cell suspension over a period of 3 
weeks. The rabbits were bled on the seventh, eighth and ninth days after 
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the last injection. Its complement was inactivated by heating at 57° for 
30 minutes. 


TABLE 1 


CoMPOUNDS INVESTIGATED 


I. HiAs€ _NHCOCONH( _AsOsH, oxanilide-p,p’-diarsonic acid 


OH 

II. Hs, ; : 
2-methyl-4,6-di-(p-azophenylarsonic acid)phenol 
R 
OH 
R, 
IIT. , 1,8-dihydroxy,2,4,6-tri-(p-azophenylarsonic acid)benzene 
OH 


R 
IV. C > AOutt phenylarsonic acid 
R represents —NNC AsO: 


Complement serum was pooled from the blood of four guinea pigs. 

Complement fixation tests along with various controls were carried out 
by mixing saline solution, guinea pig complement (1/10 dilution) and hap- 
tene or antigen solution (20 yug./ml., pH 7.0) as shown in table 2. The 
mixtures were incubated for 30 minutes at 37° and the precipitates were 


TABLE 2 
PROTOCOL FOR COMPLEMENT FIXATION TESTS 


ANTIAZOBENZENE 
COMPLEMENT -ARSONIC 
1/10 


ACID ANTIGEN OR HAPTENE 
TEST DILUTION SERUM SALINE (20 c./mc.) 
SOLUTION ML. ML. ML. ML. 
1 1.5 0 3 0 Control on complement 
2 1.5 1.5 1.5 0 Control on anticomplemen- 
tary action of antiserum 

3 1.5 0 1.5 1.5 1 
4 1.5 0 1.5 1.5 II | Control on anticom- 
5 1.5 0 1.5 1.5 III plementary action 
6 1.5 0 1.5 1.5 iV of substances I to IV 
i 1.5 1.5 0 1.51 
8 1.5 1.5 0 1.5 II | Complement fixation 
9 1.5 1.5 0 1.5 III tests of substances 
10 1.5 1.5 0 1.5 IV Ito IV 
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removed from tubes 7,8 and 9. To 0.5 ml. portions of each solution there 
were added 0.5 ml. portions of mixtures of equal volumes of 5% sheep cells 
and amboceptor of dilutions 1/25, 1/50, 1/100 and 1/200. The tubes were 
then incubated for 15 minutes at 37°. The results of the tests are shown in 
table 3. 


TABLE 3 
RESULTS OF COMPLEMENT FIXATION TESTS 


Concentration of Amboceptor 


TEST 


SOLUTION 1/25 1/50 1/100 1/200 
1 44 4+ + - 
2 4+ 3+ + - 
3 44+ 3+ + - 
4 at 3+ - 2+ + 
5 44 3+ 2+ + 
6 4+ 3+ 2+ + 
7 2+ - - 

8 2+ - - 
9 2+ - - - 
10 4+ 3+ + - 


4 + means complete lysis 
— means no lysis 


Discussion and Summary.—The experimental results show that comple- 
ment is removed from solution by the formation of a precipitate by a 
multihaptenic simple substance and its haptene-homologous antiserum. 
A simple univalent haptene, phenylarsonic acid, was found not to fix 
complement under the conditions investigated. 

It is accordingly not necessary for fixation of complement that the 
antigen be a protein or other very complex molecule; but it is indicated 
that combination with antibody (as by haptene) without precipitation 
may not be sufficient for complement fixation. 


1 (a) Landsteiner, K., and van der Scheer, J., Proc. Soc. Exp. Biol. and Med., 29, 747 
(1982); Jour. Exp. Med., 56, 399 (1982); (b) Pauling, L., Campbell, D. H., and Press- 
man, D., Proc. Nat. Acad. Sci., 27, 125 (1941); (c) Pauling, L., Pressman, D., Campbell, 
D. H., Ikeda, C., and Ikawa, M., Jour. Am. Chem. Soc., to be published. 
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MAXIMAL SYLOW SUBGROUPS OF A GIVEN GROUP 
By G. A. MILLER 


DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 


Communicated January 19, 1942 


Suppose that a given solvable group G contains maximal Sylow subgroups 
of at least three different orders. Since G is solvable it contains an invari- 
ant subgroup of prime index p. This subgroup involves all the Sylow sub- 
groups of G except those whose orders are powers of p. As the former Sy- 
low subgroups are conjugate under this invariant subgroup as well as under 
G each of them is transformed into itself under G by p times as many opera- 
tors as under this invariant subgroup. Hence none of these Sylow sub- 
groups can be maximal under G unless it is of index p under G. That is, 
this invariant subgroup cannot contain maximal Sylow subgroups of G of 
more than one order. Hence a solvable group cannot contain maximal 
Sylow subgroups of more than two different orders and if it contains maximal 
Sylow subgroups of two different orders at least one of them is invariant and of 
prime index under the group. 

If a group G contains maximal Sylow subgroups of exactly two different 
orders and two of them are invariant then they have only the identity in 
common and G is their direct product. That is, G is then the cyclic group 
of order pg, p and qg being distinct prime numbers. If only one of these 
Sylow subgroups is invariant under G then it must be abelian since other- 
wise its commutator subgroup would be transformed into itself by an opera- 
tor which would be prime to its order and G would involve maximal sub- 
groups of more than two different orders according to the theorem at the 
close of the preceding paragraph. For similar reasons the given invariant 
maximal subgroup must then be of type 1” and must admit an auto- 
morphism of prime order which does not transform any of its proper sub- 
groups into itself. This condition is obviously sufficient as well as neces- 
sary. 

From what precedes it results that a solvable group contains either no 
maximal Sylow subgroup or all its maximal Sylow subgroups are either of 
the same order or of two different orders. The last of these three cases is 
the simplest and was considered above. When all the maximal Sylow sub- 
groups of a solvable group are of the same order and the group contains only 
one such subgroup this subgroup must be invariant and hence it must be of 
prime index under the group. It must therefore involve an operator whose 
order is equal to this index and transforms into itself at least one of the 
proper subgroups of the given invariant subgroup since otherwise the group 
would contain maximal Sylow subgroups of two different orders. The 
two abelian groups of order 12 are illustrations of solvable groups in which 
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all the maximal Sylow subgroups are of the same order and invariant, 
while in the dihedral and the dicyclic groups of this order the maximal 
Sylow subgroups are of the same order but non-invariant. 

When every maximal subgroup of a given group, which may not be 
solvable, is a Sylow subgroup of G and G contains only two maximal sub- 
groups then each of these subgroups is invariant under G and G is their 
direct product since they could not have any common operator besides the 
identity. Moreover, each of these two maximal subgroups is of prime 
order since otherwise G would involve a maximal subgroup which would 
not be a Sylow subgroup. That is, when every maximal subgroup of Gisa 
Sylow subgroup and G contains only two maximal subgroups then G is the 
cyclic group of order pg, p and q being distinct prime numbers. When 
every maximal subgroup of G is a Sylow subgroup and G contains more 
than two maximal subgroups then it is not possible for all of these sub- 
groups to be invariant under G since otherwise G would contain a maximal 
subgroup which would not be a Sylow subgroup. Hence at least one of 
these maximal subgroups appears in a complete set of conjugate subgroups 
of G involving more than one subgroup, and since all the Sylow subgroups 
of the same order contained in G appear in the same set of conjugates under 
G there is one and only one set of conjugate subgroups of G which involves 
this particular Sylow subgroup. 

The given set of m conjugate subgroups is transformed under G according 
to a primitive permutation group since the subgroup composed of all the 
permutations of this group which omit a given letter is maximal. This 
primitive permutation group satisfies the condition that each of its sub- 
groups composed of all its permutations which omit a given letter is of 
degree n — 1 and that all of its permutations besides the identity are of 
degree n — 1. For if two such subgroups had more than the identity in 
common then these common operators would be transformed into them- 
selves by operators of each of these subgroups which are not common to the 
two subgroups since a proper subgroup of a prime power group is always 
transformed into itself by operators of this group which do not appear in 
this subgroup. It therefore would result that G would contain a maximal 
subgroup which would not be a Sylow subgroup since the proper subgroup 
composed of all the operators of G which would transform into themselves 
these common operators of the two given subgroups would either be a 
maximal subgroup of G or it could be extended by operators of G so as to 
obtain a maximal subgroup of G. 

As this is impossible the given primitive group of degree is of class m — 1 
and hence it contains an invariant regular subgroup of order m according to 
a well-known theorem due to G. Frobenius.! This primitive permutation 
group is simply isomorphic with G since all the maximal subgroups of G are 
supposed to be Sylow subgroups of G. For the same reason the given in- 
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variant subgroup of order m is abelian and has an order which is a power of a 
prime number. It is of type 1” since otherwise it would contain a charac- 
teristic proper subgroup and hence G would contain a maximal subgroup 
which would not be a Sylow subgroup. This proves the following theorem: 
If every maximal subgroup of the group G 1s a Sylow subgroup then the order of 
G is divisible by two and only two distinct prime numbers and there is at least 
one invariant maximal subgroup in G. If G contains more than one invariant 
maximal subgroup it contains exactly two such subgroups and is the cyclic 
group of order pq, p and q being distinct prime numbers. If G contains only 
one invariant maximal subgroup this is abelian, of type 1”, and of prime index 
under G. 

From this theorem it follows that when every maxtmal subgroup of a 
given group is a Sylow subgroup thereof then each of its Sylow subgroups 
is also a maximal subgroup. We proceed to consider the case when every 
Sylow subgroup of a given group is a maximal subgroup thereof. If sucha 
group contains a maximal subgroup which is not also a Sylow subgroup this 
maximal subgroup cannot contain a Sylow subgroup of the group since this 
Sylow subgroup is supposed to be maximal. Its index under the group 
would therefore be divisible by each of the prime numbers which divide the 
order of the group. In particular, it could not be invariant under the 
group since a maximal invariant subgroup of a group is of prime index 
under the group. That is, if every Sylow subgroup of a group is a maximal 
subgroup of the group then the group cannot contain an invariant maximal 
subgroup unless it is also a Sylow subgroup, and if it is a Sylow subgroup 
it is of prime index under the group and comes under the theorem noted in 
the preceding paragraph. 

Suppose that the group G contains maximal Sylow subgroups of order p?, 
p being a prime number. If two of these subgroups have only the identity 
in common then every two of them have only the identity in common, and 
hence they are then transformed under G according to a primitive permuta- 
tion group which contains a regular invariant subgroup whose order is 
equal to the number of these subgroups of order p* and hence its order 
would be prime to p. As the corresponding quotient group would be of 
order p* it would not be possible for every Sylow subgroup of G to be maxi- 
mal. If two of the maximal subgroups of order p* would have a subgroup 
of order p in common this subgroup would be invariant under G and to a 
Sylow subgroup of the corresponding quotient group whose order would be 
prime to » there would correspond a subgroup of G which would not be a 
Sylow subgroup of G but would be either a maximal subgroup of G or could 
be extended to such a subgroup. 

Hence it results that if every Sylow subgroup of G is maximal and G 
contains a Sylow subgroup of order p” then it contains only one such Sylow 
subgroup and hence is of order p’g, p and gq being distinct prime numbers, 
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and it contains an invariant subgroup of order p?. This invariant subgroup 
is non-cyclic and none of its subgroups of order / is transformed into itself 
by an operator of order g. The tetrahedral group and a group of order 75 
which involves the non-cyclic group of order 25 are instances of such a 
group. On the other hand when the given invariant subgroup is the non- 
cyclic group of order 9 operators of orders 2 and 3 in the group of iso- 
morphisms would transform into itself a subgroup of order 3 contained in 
this invariant subgroup and hence not every Sylow subgroup would be 
maximal. It therefore results that if every Sylow subgroup of a group is 
maximal and the group contains maximal subgroups which are not Sylow sub- 
groups then the group must not only be insolvable but each of its Sylow sub- 
groups must have an order which is at least the cube of a prime number. 


1 Cf. Miller, Blichfeldt, Dickson, Finite Groups, second edition, p. 113 (1938). 


AN INVERSE PROBLEM CONCERNING A CHAIN PROCESS 
By I. OPATOWSKI 
INSTITUTE OF TECHNOLOGY, UNIVERSITY OF MINNESOTA 
Communicated January 19, 1942 


The processes in which one or more entities change their characteristics 
in a successive manner are chain processes. They have been studied in a 
deterministic and in a statistic manner. 

In the deterministic treatment the entities involved in the process are 
grouped into a number of different kinds: {kind 0, kind 1, ..., kind }. 
Let N;,(t) be the amount of the kind 7 at the time?¢. If at the moment ¢ = 
0, in which the process started, all entities were of kind 0, then: 


Y.(0) = 1, Y,(0) = 0 fori 2 1, where Y;,(¢#) = N,(t)/No(0). (1) 


If each kind i changes directly only into the kind 7 + 1 and if the amount 
of the kind 7 changed during dt into the kind 7 + 1 is proportional to N;(¢), 
then: 


dY;/dt = k,Y;~,(t) — kp+1Y;,(t), @ = 0, 1,..., ), (2) 
where k;’s are constants with kj) = 0. If r, constants k;’s are equal to 
each other, the symbol K, will be used instead of k; (9 = 1,..., m; Dele = 
n+ 1). 

The system (2) as a mathematical expression of a chain process was 
introduced by Rutherford! to study the radioactive transformations, and 
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has been applied since then in many other fields. In chemical kinetics it 
represents a general type of successive reactions.” In the theory of bio- 
logical effects of radiations it gives a relation between the number of killed 
microérganisms and the time of irradiation.* In the theory of photographic 
emulsions‘ it explains the formation of a latent image in relation to the 
time of exposure. 

In the statistic treatment the process is a particular Markoff chain. All 
entities involved in the process are considered as forming a system which 
may take different states: {state 0, state 1,..., state m}. Let ,;(t) be 
the probability of being the system at the time 7 + ¢ in the state J, if it is 
at the time T in the state 7. If the process is homogeneous with respect to 
the time, i.e., if p,;(t) is independent of T and if ;;(t) is continuous on the 
right at ¢ = 0, then by a theorem of W. Doeblin® 9;;(¢) is necessarily a dif- 
ferentiable function when ¢ > 0 and is differentiable on the right at ¢ = 0. 
We consider a process of this type and restrict ourselves to the case in 
which an arbitrary state i — 1 may change during dt only into the state 1, 
and the probability for this is k,dt, where k; (an intensity® of the process) 
isaconstant. Then’®® (1) and (2) are satisfied by po ;(t) = Y;(#) and also 
by the other probabilities ;;(t). 

This statistic form of the process was first considered by H. Bateman® 
in connection with bombardment of a fluorescent screen by a particles. 
It occurs in general in radioactive disintegrations and was applied by W. 
Feller! to the growth of biological populations. Into the same form may 
be included also a chain process whose intensities depend on the time and 
are functions of the type: &;(t) = f(z)y¢(t), because it reduces to a process 
with constant intensities by using a suitable time variable (of such a type is 
for instance the Pélya stochastic process‘). 

In many applications the only kinds or states of the changing system 
accessible to quantitative measurements are the final ones. For instance 
in killing microérganisms by irradiation, the various degrees of injury 
which the organisms undergo before dying, are usually not distinguishable 
with an accuracy necessary for a numerical counting, whereas dead organ- 
isms are clearly identified. The same situation occurs in radioactive 
transformations when the intermediate products are unstable or when their 
radiation energy is too small for measurements.'! In such cases Y,(t) is 
the only Y;-function which is obtained from the experiments and the prob- 
lem arises to determine the laws governing the process, i.e., to evaluate 
the constants k;’s and the number ” knowing Y,(¢). This problem is the 
object of the present paper.’ Essential tools for its solutions are the 
Laplace transformation, which Bateman‘ used to solve (2), and the com- 
plete homogeneous symmetric functions‘ defined by: ho = 1; h,(Ri,... 
kn +1) = Ti = tt 224, where > is taken for all sets of non-negative integers 
a; such that }$="%*+1a; =» 2 1. Since this does not involve any compli- 
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cation we will take k,’s as complex numbers. We put also k; ¥ 0 for 1 < 
i < n, because’ if a k; = 0, all Y; =O forj 2 i. For |z| < |g>"|, the expan- 
sion :\4 


We =i(l — gz) = + Bg oy(Golp; p = 1,..., m)z, (3) 


holds, where h,(q,|7,; p = 1,..., m) means that each q, is taken r, times in 
calculating h,. We need an expansion in nen fractions: 


Bot(s + K,)~* = Destro 'C,(s+K,)*~. ©) 


Introducing here s = z — K,, K, — K, = q,° for p ¥ 1, multiplying by 2” 
and taking into account (3) we get: Doi”, ~-(—q,)(1 — g,2) * = 
eso" 'C,,2" + terms in z of degree 2 r,. Therefore, by (3): 

C., = hot a(K, — K,)” 
where /, stands for 

h,((K, — } Se hate” a= 1,.5.)9j. 9 9). 
The Laplace transform" of Y,,(2), i.e., y,(s) =@{ Y,,(t), s} =U Y,} is: 
yn(s) = Soe” “Y,(t)dt = PIGS i(s + K,)~’, 
where P,, = kike...Ry. (5) 
Therefore by (4), since! v!(s + K)~"~' = Q{fe~**'}: 
V,{0) = PSS Bese gfe 7, 
where A,, = C,,/(r, —u — 1)! (6) 


If ally, = 1, then A,, = A,,9 = C,, = C,,o = Witt .(k — &,)7', by 
a formula of Bateman." It is seen from (6) that!® 


Q{ Y,,(t), s} exists for Rs > —RK,, p = 1,..., m. (7) 


Y,,(t) may be written for nm 2 1 asa repeated convolution: 
Y,(t) = P, exp(—Ait)*xexp(—fot)«...xexp(—R, + 14), (8) 


which is equivalent to a repeated integral, because e~™s F(t) = e~™ fo'e™ 
F(t)dt. It is seen from (6) or (8) that Y,,(¢) is a symmetric function of ki,..., 
k,. Putting in (3) q, = —K,,2 = 1/s, we get: s” + 1 = (s+XK,)~* = 

Alt o(—1)"h,(K, |r, ; p = 1,..., m)s~” for |s| > |K,|. Therefore by (5): 
yn(S) = Padi = of — iy tte ae?  ° ~), whee hlhy..c hed 
=h aN Ir; p = 1,..., m), and by known theorems," since 2{f*} = a! 


-a- 


s 
Va(t) = Padir'= (—1)"y(hi,---s kn dE *"/(v +m) (9) 
which" by (6) or (8) is convergent for any finite ¢ and ;. 
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We give now three methods for the solution of our problem: the deter- 
mination of the k,’s and of the m from Y,(¢). The question what method 
is more convenient in a practical application depends on how Y,(t) is 
known. If it is known with a very good accuracy for small values of t 
method I may be chosen. Method III requires an accurate knowledge of 
Y,,() for large values of ¢ if n is large. 

Method I. We have from (9): d*Y,(0)/dt® = Y“(0) = Oforx <n, 
yi" +0) = (—1)’P,h,(hi,..., Re 41) for vy 2 0. Since ho = 1, is ob- 
tained as the order of the first derivative of Y,(¢) which is ¥ 0 at t = 0. 
Besides this: 


Iey(Iess.. 2» Be +1) = (—1)° YE * 9 (0)/¥L"(0). (10) 
Therefore hi,...,k, 4 , are determined as the roots of the equation 
x8 a +... ta,41= 0, (11) 
where (— 1)'a; are the elementary symmetric functions of hi,..., ky 4 , and 


may be calculated by means of the equations )% = fh; — yy = 0, a0 = In 
= 1, 
Method II. This method and the following one hold if Rk; > 0 
fori < nand k, 4 ; = 0, orif 
Rk; > Ofort KC n+ 1. (12) 
(These restrictions on k,’s are normally satisfied in all applications.) We 
assume first that (12) holds. Then by (6): Y{(o) = Oforx« > Oand 
So? VS (Odt = —Y&- (0) for x 2 1. (13) 
Therefore n + 1 is the smallest value of x > 0, for which the integral (13) 
is #0, and by (10): 
h,(ley..-1 be +) = (1K VE Tt (at/ VE * dt. (14) 


When h,’s are known, &;’s are calculated as in the method I. It is possible 
to lower the degree of equation (11) by one. In fact, by (7), (12) and (5): 


S{V,(, 0} = yn(0) = So? Valt)dt = 1/ky +1, (15) 
which gives k, +1. Then" 
h,(Ra,..+5 ky) = h,(hi,. oe) ky + 1) aaa ky + ih, - i(i,.. ” Rn ~ ») 


may be calculated by means of (14) or (10) for vy = 1,..., m, and hy,..., Ry 
are determined as the roots of an equation of type (11) of degree n. 

As soon as k, 4 ; is known from (15), 2 may be determined also from the 
formula 


n = 1+ lim,-, offd log.|¥,() + kn+1¥n(#)|/di], — (16) 
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which holds without any restriction on k,;’s. To prove (16) write a power 
series for Y,, — ,(t) by changing in (9) m intom — 1. From this it is easy to 
see that 


lim, _, olt Vy ~1()/Y, -1@) =n" — 1. (17) 


On the other hand (2) gives forz = n: log| Y;,(t) + Rn + iY, (¢)| = log ln + 
log LY, “ i()|. By differentiating this with respect to ¢ and taking into 
account (17) we get (16). 

We show now that the method remains substantially unchanged if k, + ; 
= (and if (12) holds only fori < n. Ofcourse (16) is true and all the other 
formulas hold too if m is changed in them into  — 1. But since by (2) 
Y(t) = RknY, — (8), (13) holds also for Y“ if x > 2. (15) gives now 
So? Vn —i(t)dt = 1/k,. Therefore (t) Jo° V;,(t)dt = Rn Jo” Yn — 1(t)dt = 1 
and » + 1is the smallest value of x 2 2 for which the integral (13) is ¥ 0. 
XV; (t) = 0, by (2), because ky + , = 0, therefore by (1) }Y,(t) = 1. On 
the other hand from (f) and (1) we get Y,(©) — Y,(0) = Y,() = 1 for 
n > 0, and we conclude that the system which was entirely in the state 0 
att = Oisatt = o entirely in the final state n. 

Method III (Method of Moments). We assume first that (12) holds. 
Then by known theorems of the Laplace transformation :" 

ye(0) = (-1" RY, 0} = (-V’K°VWde (18) 
is convergent. On the other hand from (5): Ry 4 wWa(s) = ms?t i+ 
k;-'s)~* which by (8) is = °- o(—1)"h,(ki',..., ka + 08” for |s| < [ez"]. 
Therefore: ky + (0) = (—1)’vlh,(k7',..., Ra 4. 1) and by (18): 


rh (h;',..-. feed = bnarthe FYOM, » = 0,1,..: (19) 


This reduces the determination of the k;’s to the solution of an algebraic 
equation in the same way as in the method II. 

If (12) holds only for i < n and k, + ; = 0, all formulae of this method 
hold if we change in them into n — 1. Therefore, since now Y,(t) = 
knYy —(t), (19) becomes: h,(Rk7’,..., ka’) = ("Wo °rY, (dt, which 
determines ;,..., ky, in the same way as previously.” 


1 Rutherford, E., Chadwick, J., Ellis, C. D., Radiations from Radioactive Substances , 
Cambridge, 1930, pp. 11-12. 

2 See for instance Branch, G. E. K., Calvin, M., Theory of Organic Chemistry, New 
York, 1941, p. 357. 

3 Crowther, I. A., Proc. Roy. Soc. London, B100, 401 (1926); Opatowski, I., Bull. 
Amer. Math. Soc., 47,704 (1941); 48, 46 (1942). 

4 Cf. for instance Silberstein, L., Jour. Optical Soc. 31, 343 (1941). 

5 Doeblin, W., Bull. d. Sc. Math., 62, 21-32 (1938); 64, 35-37 (1940). 

6 Lundberg, O., On Random Processes and Their A pplication to Sickness and Accident 
Statistics, Uppsala, 1940, pp. 57-68. 

7 Kolmogoroff, A., Math. Ann., 104, 428-437 (1941). 
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11 Meyer, S., Schweidler, E., Radioaktivitaet, 2nd ed., Berlin-Leipzig, 1927, p. 32. 
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sum of exponentials Dc; exp (—k,t) = Y,,(t). In the case in which ¢;’s are independent of 
k,’s various fitting methods are known: Aigner, F., Flamm, L., Phys. Ztschr., 13, 1151- 
1155 (1912); Walsh, J. W. T., Proc. Phys. Soc. London, 32, 26-30 (1919); Levy, H., 
Ibid., 34, 108-113 (1922). Further references in Meyer-Schweidler! p. 61; also 
Bernstein, F., Ztschr. ang. Math. u. Mech., 7, 441-444 (1927) and Doetsch® p. 27-28. 
In our problem however c¢,’s are certain functions of k,’s, and these methods are not good. 

13 Bateman, H., Proc. Phil. Soc. Cambridge, 15, 423-427 (1910). This paper remained 
unknown to some recent authors who gave less simple methods for the solution of (2): 
Arley, N., Mat. Tidsskr., 49-51 (1939); Lundberg;® Feller.’ 

14 Murnaghan, F. D., Theory of Group Representations, Baltimore, 1938,pp. 108, 112. 
These functions are called also homogeneous product sums; cf. Littlewood, D. E., Theory 
of Group Characters, Oxford, 1940, pp. 82, 83, 88 or, alef-functions of Wronski, s. Enc. d. 
math. Wiss., IB3b, 465, 459. 

18 Doetsch, G., Theorie der Laplace-Transformation, Berlin, 1937, pp. 23, 148, 157, 43- 
48, 61-64 and Widder, D. V., The Laplace Transform, Princeton, 1941. 

16 Only the first term of this expansion was known, s. Handb, d. Phys., 22, 190. 

17 The author acknowledges helpful information received in preparing this paper from 
Profs. J. L. Doob, L. Silberstein, J. W. Tukey. 


ON CONFIDENCE INTERVALS 
By EpwIn B. WILSON 
HARVARD SCHOOL OF PUBLIC HEALTH 
Communicated December 23, 1941 


In 1927 I called attention to the fact that many statements about proba- 
bility are highly elliptical and illustrated the matter by the simple case of a 
point-binomial universe with unknown probability p and observed value po 
in some sample.! Using the admittedly rough estimate of probability 
based on the standard deviation one ordinarily writes 


bo — > V pogo/n < p < po + 4 V pogo/n 


and states that the probability that the true value » in the universe lies 
between the limits given may be had from a probability-integral table 
entered with a normal deviation of A\units. I urged that a better procedure 
would be to use for the standard deviation the value ~/ bq/n obtained from 
the unknown of the universe which leads to 


era yr erg 
po+t/2 Vogt +H/4  , . oti/2 poate + F/4 1) 








1+? i+ 1+ 1+ 
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where ¢ = \?/n, and to the statement that if (unknown) should be outside 
of those limits the probability of observing the sample would be less than 
that obtained from entering the probability-integral table with a normal 
deviation of \ units. I was trying to emphasize that we know nothing 
about the value of », which must have whatever value it did have in the 
universe from which the sample was drawn, but that we could set limits 
based on probability calculations such that if p lay between them the chance 
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Confidence levels in terms of 1 — P = 1 — 2e for n = 10, x =9, po = 0.9 as functions 
of the unknown probability » on two assumptions: 

1 (a), left: The no less extreme observation (Clopper and Pearson), 

2 (b), right: The no more probable observation (Wilson). 


of getting the particular observation or any less probable one would exceed 
some preassigned value P whereas if p lay outside them the chance of get- 
ting the observation or any less probable one would be less than P. 

In 1934 Clopper and Pearson? introduced confidence intervals for what I 
take to be nearly the same notion that I was discussing. In their termi- 
nology, if 2e be any positive number on which to base a certain level of con- 
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fidence in terms of probabilities, the interval (1) would have a confidence 
1 — 2e = 1 — P as judged by the test I used. Their test was different; 
they asked the question: For what value p; of p less than » will the value 
po or a greater value arise with a probability « and for what value pe of p 
greater than fo will the value of p) or a smaller value arise with the same 
probability «? It may be noted that they divided the probability P = 2 
into two equal parts and asked two questions, as though one were thinking 
of deviations from a central tendency, whereas I had kept to a total proba- 
bility P and had asked a single question in terms of samples no more prob- 
able than the one observed.* Furthermore Clopper and Pearson (and 
Rietz) introduce the notion of confidence belts, which I did not have in 
mind, and state: ‘We cannot therefore say that for any specified value of 
x( = npo) the probability that the confidence interval will include p is 0.95 
ormore. The probability must be associated with the whole belt, that is to 
say with the result of the continued application of a method of procedure to 
all values of x met with in our statistical experience’’—a statement I 
should hesitate to make. 

In figure 1 will be found for n = 10 the graph of the values of 1 — 2e = 
1 — P against those of p for the case p) = 0.9 or x = 9 for (a) the Clopper- 
Pearson assumption that ¢ is the probability of all values 9 or less or of all 
values 9 or greater for x according as p > 0.9 or p < 0.9 and for (b) my own 
assumption that P is the probability of the sample of 9 or of any sample no 
more probable than it. The figure shows that on the Clopper-Pearson 
basis the curve is a continuous function of p except at the value p = pp) = 
0.9; it also shows that the shorter confidence intervals around » = 0.9 have 
confidence values (probabilities?) less than zero. The figure shows also 
that on the basis of my definition involving the probability of all no more 
probable cases the curve of P has multiple discontinuities. 

The question of generalization may be raised.* Let there be given two 
independent samples of and n’ from universes supposed to have the same 
unknown value p but furnishing from observation two values p» and py’ in 
the two samples. If the numbers in the samples are large enough to justify 
estimates of probability based on values of x?, one may write 





ne (ho — mp)? (n'po' me i! 238 (2) 
npq n pq 


x°(p — p*) = npr? + n'po'? — 2(mpo + n'po')p + (nm + n’)p?, 


which is a quadratic equation for the unknown value of p. Let 


n+n'=N, np + n'po'? = Np? + Noy’, t = x*/N; 
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then 
Pot t/2  V poget + 2/4 — of°(1 +4) po + t/2 
1+ 1+ ita 2% ba 
V pogot + 12/4 — 0,7(1 + 2) 3) 


1+ 


This result is entirely similar to (1) except for the addition of the term 
—o,?-(1 + #) under the radical; it is, furthermore, entirely general for any 
number of samples drawn from different universes having the same , 
provided $ be the weighted mean p for all samples and o,’ the weighted 
variance of the observed values fo. 

The new term —o,?(1 + #) under the radical, being negative, shows 
that for a given value of t = x?/N the range for p is the more restricted, the 
greater the scatter of the values ) observed in the different samples; indeed 
we must have® 
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Confidence levels in terms of P = 2e as functions of the unknown probability p on 
the hypothesis of the no more probable observation for two cases: 

2 (a), left: nm =n’ = 5000, po = 0.49, po’ = 0.51. 

2 (b), right: n =n’ = 3, P= 1, fo’ = 0. 


to have any range for p. If we have k samples and enter a x? table in the 
row for k degrees of freedom under any value of P to find x? and hence ¢ = 
x?/N, N being the sum of the numbers in the k samples, we may compute 
from (3) the limits for p such that if p lies outside those limits the chance of 
getting the set of samples and all less probable sets is less than P; if p lies 
within the limits (2) the chance is greater than P but is in any case less 
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than the value of P corresponding to that value of ¢ or x? which makes the 
radical in (3) vanish or corresponds to the equality sign in (4) and thus 
shrinks the range of p toa point.’ Figure 2a gives a plot of P against p for 
the case n = n’ = 5000, N = 10,000, po = 0.49, po’ = 0.51, po = 0.50, 
a,* = 0.0001. 

When the numbers are small the question of the test that replaces x? 
has to be discussed. If we use the criterion of the total probability of the 
samples observed and of all samples of no greater probability, we may with 
patience discuss any particular case. For example, ifm = 3, p) = 1, n’ = 3, 
po’ = 0, the joint chance of the observed samples is p*g* whatever the un- 
known value of p. There are (m + 1)( n’ + 1) = 16 possible pairs of sam- 
ples. The graph of the probability P as a function of p isin figure 2). The 
graph is, of course, discontinuous. It shows that there is no value of p for 
which the probability P exceeds 0.0625, and this for only the single value 
p = 1/2 which on a priori grounds would be infinitely rare. Values other 
than p = '/2 between 0.37 and 0.63 give no values greater than P = 0.047. 
There are very short ranges for p near 0.36 and 0.64 which give P > 0.05. 
Should we desire confidence intervals* corresponding to P = 0.05 we 
should have to specify something like 0.360 to 0.366, 0.5, 0.634 to 0.640; 
should we desire confidence intervals corresponding to P = 0.01 we should 
have something like 0.17 < p < 0.83; for P = 0.03 we should have 0.248 to 
0.250, 0.302 to 0.366, 0.385 to 0.615, 0.634 to 0.698, 0.750 to 0.752. The 
exact method of R. A. Fisher, if he would apply it to a problem of this type 
(as is unlikely, because the approaches seem to belong to different universes 
of discourse) and if in this symmetrical case he would consider together both 
ends of his series, would give P = 0.10; there is no value of the unknown p 
which gives a probability as much as two-thirds of this amount. It thus 
appears that the exact method affords a different criterion from that of the 
total probability of the no more probable experience, just as both of these 
methods furnish different criteria from that based on the distribution of the 
difference po — po’. 


1 Wilson, E. B., ‘“‘Probable Inference, the Law of Succession, and Statistical Infer- 
ence,’ J. Amer. Statistical Association, 22, 209-212 (1927). 

2 Clopper, C. J., and Pearson, E. S., ‘““The Use of Confidence or Fiducial Limits Illus- 
trated in the Case of the Binomial,’’ Biometrika, 26, 404-413 (1934). An account of 
this may be found in Rietz, H. L., ‘A Recent Advance in Statistical Inference,’’ Amer. 
Math. Monthly, 45, 149-158 (1938). 

3 There is a third notion not very different from mine or that of Clopper and Pearson— 
or rather of Neyman who slightly antedates them, at least with respect to continuous 
distributions—with which Fisher, R. A., on “Inverse Probability,’’ Proc. Camb. Phil. 
Soc., 26, 528-533 (1930), associated the term fiducial limits and according to him limited 
to continuous distributions, see ‘‘The Fiducial Argument in Statistical Inference,’ Ann. 
Eugenics, 6, 391-398 (1935). 

4 If I interpret Clopper and Pearson aright, one may logically take 2e = 1, e = 0.5 
and construct a confidence belt, quite as one constructs one for 2e = 0.05. The belt is 
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of course thinner than for smaller values of ¢ but, still, has an appreciable area within 
which the probability that p will lie as ‘‘the result of the continued application of a 
method of procedure to all values of x met with in our statistical experience”’ is zero. 
As I do not see how this can well be, I incline to some doubt in respect to their (or Rietz’s) 
statement. Moreover, for my test as shown in figure 1 (5) there is a finite interval for 
which P = 1 and the confidence is 0, whereas the probability that p lie in that interval 
could hardly be 0. Thus although confidence intervals are based on probabilities, it is 
not certain that probabilities are based on them. The difficulty of the negative values of 
probabilities 1 — 2¢ for p lying in intervals close around pp = 0.9 can be obviated by a 
modification of the definition whereby when p > fo we use only a fraction @ of the maxi- 
mum term M for x = nmpyand when p < fo the fraction (1 — 0) of M writing X + 0M =6e 
and (1 — 6)M + Y = ewhere X is the sum of the terms up to but not including M, and Y 
is the sum of the terms beyond M, the value of @ being !/, + 1/.(Y — X)/M. 

5 [ shall not try to offer methods of generalizing the Clopper-Pearson test. 

6 While o,,? may in extreme cases be nearly equal to Pog so that x2 may be near 2No, = 
2N(Pogo)'’?, the interesting range for o,? when N is large is a small multiple of pogo/ N 
and thus ¢,? is small compared with Pogo, and approximately 


ee. i PN. oe Oper 
srikse ogo/N Pogo 4 pogo} | 


The expression in front of the bracket is the usual expression for x? taken from the 
marginal totals of a 2 by k table and the bracket is nearly 1. Throughout this paper as 
in that of 1927 I have carried out formal algebraic operations as though they were 
significant; it should, however, be borne in mind that the expression (2) for x? is dis- 
continuous and gives from a x? table the probability of the given sample or set of samples 
and of all no more probable ones only to a certain degree of approximation even when 
the numbers are large—and thus the correction in the bracket may be illusory. 

7 It follows that in the case now under discussion the probability P which corresponds 
to a certain range of p such that if p be outside it the probability of so aberrant samples is 
less than P does not correspond to a confidence 1 — P that p lie within the interval. 
If Po be the probability which shrinks the interval to a point, we might perhaps say that 
Po — P is the level of confidence that lie within the interval and 1 — Po the level of 
confidence that there be no value of common to the two universes. 

8 As in the previous case, the chance if p lies outside the interval is less than the 
specified amount, but the chance if p lies within the interval is between the specified 
value and some upper limit Po, leaving 1 — Po as the level of confidence that there be no 
value of common to the two samples. It may be repeated that the chances involved 
are not those that there actually is a common probability lying within or without certain 
limits but the chances that the observed samples and all no more probable should arise 
from such a hypothetical common value of p. 
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ON CONTINGENCY TABLES 


By EpwIn B. WILSON 
HARVARD SCHOOL OF PUBLIC HEALTH 
Communicated January 8, 1942 
In 1934, Yates,' acting on a suggestion from Fisher, developed the so- 
called exact method of treating 2 X 2 tables and proposed and discussed 
what has become known as the Yates’ correction for the x? expression used 
on such tables. Soon thereafter the method was inserted in a new edition 


of Fisher’s Statistical Methods for Research Workers (Art. 21.02) and be- 
came rapidly and widely adopted. The method consists in considering the 


linear series of c + d + 1 tables 

















a—djb+d a—1/b+1 a|b 

etd: 6 °°" e4112=-1 cl€ 
a+1\db-1 a+clb—c 
c-—l1jid+1'°°” 0 Ic+d 





where the convention has been adopted that c + d is the smallest marginal 
total, and of computing the relative probabilities 
ee ae ere Fe) 1 

(a+b+c+d)! (a — d)b+d)"c +d)! "” 











1 1 
albicid! "’ (a + c)'(b — c)!0Ne + a 
of those c + d + 1 tables. The properties of this series of probabilities, 


whose sum is 1, are therefore of great importance. Connected with them 
is the algebraic identity 


» Me Bee (a+b+c+d)! 
albicld! (a + b)"a + 0)\(b + d)\c + d)! 


where the summation extends over all the c + d + 1 terms of the series. 
The series is symmetrical when and only when two of the marginal totals 
on the same side of the table are equal.’ 

The mean value of a will be 


esse sa a 
(a+b+c+d)! 
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The question is as to the value of the second sum }>’ which is extended 
over all tables or terms. For those tables in which a = 0 the sum )> van- 
ishes as does the sum })’ because (—1)! is infinite. Hence in evaluating 
either sum one may reject the tables for which a = 0. Consider now the 


. , sie a 3 , 
series of tables arising from ——— : which will have the same values of 
‘s 


(a — 1)!d!cld. as the given series. Then 





7. oe SB oes 5 oe) 5 a ae ae ee 

(a—1)blcld! (a+b—1)"Kat+c— 1I)b+d)\c +d)! 
(a + b)(a +c) 
a+b+c+d 


where V = a + 6+ c¢ +d and where p; = (a + 0)/N, pe = (a + c)/N are 
the proportions observed in the margins. The proof may be extended to 
cover the average of any (partial) factorial as: 


a(a—1)..a@—p+1)= 


(atb+c+d—p)\a+b)\(a +)\(b +4) +a)! 
atb+c+d)at+b— p)iatec— p)b+d)\c +d)! 





and hence d = 


= Np;p», 








From the values of the factorials the moments y; of the distribution of a 
may be obtained with routine algebra. Let 


x=a— G= (ad — be)/N. 











Then? 
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ba(x) = V-1 Pipeng = Spal Noe) 
N3 
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If we set 
a Pa eed. ae a ee ad — be ed Nx 
oan” 2.40." © S£8Gte Beets 


the characteristics of the distribution of p — p’, ie., of the difference be- 
tween the proportions observed in the first and second rows of the tables 
when the tables are restricted to the linear series used by Yates and Fisher, 
may be obtained at once, e.g., 





we(p — p’) = N"Pibedde je 3g 
. (N+ 1)(a+b)%e+d)? N-1pm 


If the two rows had arisen by independent sampling in two universes with 
probabilities p and p’, respectively, and with numbers” = a + b,n’ =c+d 
in the two samples, there being (m + 1)(m’ + 1) different complexions 
for such pairs of samples, the value of ue for the observed values of p — p’ 
would be pg/n + p’q’/n’ and if it should be assumed that p and p’ were 
equal and had the value (a + c)/N obtained from the marginal totals we 
should get uo(p — p’) = poq2/(Npigi) which is slightly less than that ob- 
tained above. If it were assumed that the 2 X 2 table had arisen by 
sampling from a non-associated universe with probabilities pipe, pige, Pog, 
992 for the four cells occupied by a, }, c, d, respectively, the relative fre- 
quencies p, p’ in the two rows as observed could have no frequency distribu- 
tions because NV + 1 of the (VN + 3)(N + 2)(N + 1)/6 samples would turn 
out witha + 6 = Oand N + 1 of them with c + d = O and the ratios for p 
or p’ would be 0/0; the total probability that p — p’ would be indeter- 
minate would ordinarily be small but the fact that there would be indeter- 
mination is sufficient to show that there would logically be no distribution 
for p — p’ and hence no moments.‘ However, the expression x = (ad — bc)/ 
N, though no longer interpretable in terms of the difference p — p’ 
would have a distribution, not of c + d + 1 elements nor of (a + b + 1) 
(c + d + 1) but of (V + 3)(N + 2)(N + 1)/6. The value of ue for this 
distribution is readily found to be® 


po(x) = (N — 1)pipegige compared with Npipogigi: or ear Pibegige 


when the sampling is by the double point binomial or by the Yates-Fisher 
series, respectively, the values of the unknown probabilities in the uni- 
verses in question being taken from the marginal totals. 

The argument by which the linear series of tables with their relative 
probabilities are obtained from a given 2 X 2 table by considering varia- 
tions in the body of the table consistent with the given marginal totals is 
equally applicable to m X n tables. Given the table 
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if II be the sign of a product, is its relative probability in the (m — 1)(m — 1) 
dimensional spread’ of tables derived therefrom by all variations of aj; 
consistent with the marginal totals A; and B;. The mean value in such a 
spread for the term a;; is A;B;/N. The average value of a;;(a;; — 1) or of 
A;Qip, k # j, or Of a,x, k ~ 1,1 # j, can be found by an argument en- 
tirely similar to that given before. For example, 


aise i DAs > : Qy234 > ‘ 1 | 
a 2 eS i ae 
1234 N! I1;;a;;! (aw — 1)\(as4 — 1)!11’a;;! 


where ),’ extends to all variations consistent with the former marginal totals 
except that A; is replaced by A; — 1, B; by Bz, — 1, As by A; — 1, By by 
B, — 1 and II’ extends over the mn — 2 elements other than ay, and ayy. 
A:BrAsBy 

N(N — 1) 

ne A;BrA;B,  AiBe A3By 
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or the mean product relative to the mean is 
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N? é ; 
7 (AjjQp,) = way Did; Peg, k A 1,1 Aj 
where p; = A;/N, 9; = B;/N. Similarly 


N? 
be(ayj) = par pil — pg — gy) 


N2 
7(4;0;,) = — fare Pill — Pig 


Hence the correlation coefficients of two terms not in the same row or 
column and of two terms in the same row but not in the same column are 


ea PiDiPrQ 
r(ajj~) = + Vi — p)(1 — g)(1 — Pa) (1 — g)’ 


gel qi 
1(ajjQi) = Vi — qg)(1 — aq)’ 
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Higher moments could be obtained, with sufficient patience, if they were 
needed. 

It sometimes happens that the numbers in two frequency functions are 
really small. For example, the distribution by age of deaths from chicken 
pox for the two sexes as reported in Massachusetts for 1934 were 


Age <1 12 3 4 5-9 10-14 15-19 20-29 30-39 40-49 50-59 Total 
Male Ba Say op 1 0 0 0 1 0 1 7 
Female 02102 0 0 0 0 0 0 0 fs) 
Total So ae oe 0 0 0 1 0 1 12 


with none for all higher age groups. This as it stands without the higher 
age groups is a 2 X 12 table; but it seems rather clear that it has to be 
considered as a 2 X 7 table, omitting the columns with no elements.® 
If so considered x? = 9.28 looked up under 6 degrees of freedom gives P = 
0.17 which is not significant. There are 130 different tables consistent 
with the marginal totals. The relative probability of the given table is 
1/264. The 130 relative probabilities which arise have the frequency 
distribution 


1 2 3 6 9 18 

prob. = eae — a4 a ——e 
792 792 792 792 792 792 

freq. 18 9 40 22 26 15 


The total probability of the given table and all no more probable is 52/264 
= 0.20 which is not much larger than the value previously obtained from 
x?. However, the variance of the age for the males is 478 with 6 degrees of 
freedom and for the females is 2.4 with 4 degrees of freedom which gives 
a ratio of 208 and the 0.01 point is at 15.2. It is far from certain that 
variance analysis can be relied upon; but almost any test for the signifi- 
cance of the standard errors of the distributions would indicate a high 
degree thereof—thus indicating how different the result of different tests 
corresponding to different formulations of the problem may be. 


1 Yates, F., ‘Contingency Tables Involving Small Numbers and the x? Test.’’ Suppl. 
Jour. Roy. Statist. Soc., 1, 217-235 (1934). 
n—k| k n—k| : : 
2 The tables ——— | ——— and ——— | ——————- where k < '/gn have the interesting 
k n—k k |In-k—1 
property that the total probability of the given tables and of all tables in the respective 
series which are no more probable than these (taking both tails together) are equal and 


equal to 


2(nl)? f ae n(n — 1) 3 n(n — 1)...(n —~ k + 1) | 
ny! ‘gas +| yt | tot] i it 











or 
oul 
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a a?| 4, ther have th ililihins tian 1 esa the ilies 
=e 1 er have e same proba == 25 i - 2 e ta — 
mig - ¥35 013 


1 6/0 | 


—" 


. . 
3 | 0 have the same probability as : I’ ols’ 5 


0’ namely, 37/462. 
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3 As an illustrative example consider the table a = 10, b = 3,c = 2,d = 15, N = 30. 
The 18 relative probabilities areas 1; 102; 2992; 37,400; 235,620; 816,816; 1,633,632; 
1,925,352; 1,337,050; 534,820; 116,688; 12,376; 476, the common denominator being 
6,653,325. The series looks skew in the arithmetical sense, but actually « = 0.02113; 
also 6. — 3 = —0.08084. If we note that for a point binomial the skewness « is (¢ — p)/ 
/ Npgq we note that for x in the Yates-Fisher series x = ~/N — 1(1 — 2/N)xix2 and thus 
varies for large values of N jointly as the square root of N and the product of the skew- 
nesses of the point binomials in the margins. 

‘If N = 100, and p, = 0.02 for the universe, p. being anything, there would be 2 ele- 
ments on the average in the first row which would, however, be entirely vacant in e~? or 
13.5% of the samples. 

5 Had we assumed a general four-fold universe the average value of x = (ad — bc)/N 
and of its variance would have been 


z 2 
x = (N — 1)(m — pifr), ua(x) = (N — 1) Es — (: =_ = &? + 


(: = » 5(q. — fi)(Q2 — p» | 


A - ™ 
age ae nae a 
6 Let it be assumed that qi, go, ..., Ym are the probabilities that an element fall in 
column 1, 2, ..., m, respectively, in any row. Then the probability that ay, ai, ..., 
dim, respectively, will fall in the cells of the first row is 


where 


Ti | 





56 = m — pip: and the universe is 


Aj!q:°"q."" oes Imm wil 11,A;! Bi, Bs 
44; !ay9! wae Aim! Il,a j! ne o-+ Ge 


Bm 








will be the probability of the m X mn observed frequencies in the table. The unknown 

factor involving g; will be the same for all tables having the same marginal totals. But 

the probability of those marginal totals arising from the distribution of N elements is 
N! T1;A ,111;B;! 


—— q,3ig,P* ... g,,2m and hence the quotient 
11,B;! qi Qe Im q Nya! 


is the relative probability for any distribution a;; consistent with the marginal totals. 


7 Whether ‘‘dimensions”’ or ‘“‘degrees of freedom”’ or ‘“‘spread’’ are proper terms to use 
for variations when numbers are very small may be questioned. One may illustrate the 
procedure for a particular case. Consider 

3, 0, 2 5!(3!)32! 1 1 


0,3,0 and its relative probability 31 (322108 = 56 


together with all its variants and their probabilities (in brackets) 


3,1,1 [6] 3,2,0 [3] 2,12 [9 ] 2,2,1 [18]. 
0, 2, 1’ | 56 |’ 0, 1, 2’ | 56 |’ 1, 2, 0° | 56 |’ 1,1, 1° | 56)’ 


2,3,0 | 3 | 1,3,1 | 6 |] 1,2,2 | 9 | 0,3,2 [1 
1, 0, 2” | 56 |’ 2, 0, 1° | 56 |’ 2, 1, 0° | 56 |’ 3, 0,0 | 56] 
The probability of the given table and all no more probable ones is 2/56, and if one uses 


a level of 0.05 for significance, and uses the criterion of the given table and all no more 
probable ones, this table would be significant. Nothing has been said of the sampling 


























100 MATHEMATICS: H. WEYL Proc. N. A. S. 


process by which the table arose and if I understand Yates! aright it makes no difference 

whether it arose by distributing 8 elements over the six cells from a universe with the 

probabilities = | ed md with Zr = 1 or by distributing 5 elements and 3 elements, respec- 
4 5 6 

tively, to each of two trinomial universes for which m7, + 32 + 33 = m4 + 45 + ms = 1 or 

by distributing 3, 3 and 2 elements, respectively, to each of three binomial universes for 

which 7 + ™ = 2 + mr = 13 + 16 => & 

8 If the tabulation were by sufficiently fine age groups the two rows of the contingency 
table would have nothing but 0’s or 1’s and so would the totals in the columns. If 
there were 7 entries of 1 in the top row and 5 entries of 1 in the bottom row and 12 entries 
of 1 in the marginal totals along the bottom, there would be 12!/(7!5!) = 792 different 
tables all of equal probability 1/792, and the total probability of the no more probable 
tables would be 1. 


CONCERNING THE DIFFERENTIAL EQUATIONS OF SOME 
BOUNDARY LAYER PROBLEMS. II 


By HERMANN WEYL 
INSTITUTE FOR ADVANCED STUDY 
Communicated January 29, 1942 


In a previous note under the same title’ I have solved by successive 

approximations the following boundary value problem: 
‘ ( w’’’ + Qww’’ + 2(k? — w’) = O for z = 0; 

< 
(Ay) w=w’ = 0 forz = 0, w'(~) =k, 
involving the positive parameter k, for the special values \ = 0 and !/2. 
Convergence of the process for \ = '/: requires a certain constant g to be 
<1. The actual computation of q is laborious; however, I wish to report 
that a comparatively simple estimation yields g < */4. 

For arbitrary } = 0 I have to take refuge to a much more sophisticated 
method, that of fixed points of functional transformations, which was in- 
augurated by Birkhoff and Kellogg in a famous memoir? and later formu- 
lated in general abstract terms by Schauder and Leray. The experiences 
with the special cases suggest the following procedure. We eliminate the 
parameter & from the differential equation by differentiation, 


w’’’’ + 2ww’’’ + 2(1 — 2dA)w’w” = 0, 
and observe that the function x.w(xz) solves the new equation of fourth 


order if w(z) does, whatever the positive constant x. Hence we seek our 
solution of (A,) in the form w(z) = x.f(xz) where 


f+ aff’ + 21 — a)f’f” = 0 for z = 0; " 


f(0) = f'(0) = 0, f"(0) = 1; f"(o) = 0. 
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Our line of attack can be explained a little more easily if we replace the 
infinite by the finite interval 0 < z < a and thus the last boundary condi- 
tion in (1) by f’(a) = 0. We then introduce f” = g = gas the unknown 
function and start with this set-up: 

f=) = Sos — Salsas. (2) 

go” + fo’ + 21 — 2W)f’e = 0; | 
9(0) = 1, g(a) = 0. 
g = &. (4) 
For an arbitrarily given g we form (2), then solve the linear boundary 
value problem (3) of familiar type, thus defining the functional operator ®, 
which carries g into ¢, and as the last step (4), we try to determine g as a 
function whose image coincides with the original, 
g = &lg}. 


The success hinges on the following 
Lemma: Under the assumption g(z) = 0 the linear problem (3) has a 
unique solution ¢; it satisfies the conditions 


¢ 20, ¢’ + 2fe <0 (5) 


(3) 


throughout the interval 0 S z S a. 
To prove this we set 


p(2) = exp(2fo°f(s)ds) 
so that p(0) = 1, p’ = 2fp, and introduce the auxiliary function 
v1 = P(g’ + 2f¢) = (p¢9)’. 
The single differential equation for ¢ is then replaced by the system 
yg’ + fo = w/?, 
1’ — fe, = ArPf’¢. 
Multiplication by ¢, and ¢, respectively, and addition leads to the formula 


ge / 'G wt nef eas b<e). 
b 


(6) 


The inequality g = 0 which implies f’ 2 0 guarantees the positive definite 
character of the ‘‘Dirichlet integral’’ at the right side of (7). Our lemma is 
an easy consequence of this fact. 

Thus we have defined the operator , for all g 2 0; the restriction g 2 0 
is harmless because it carries over to the image yg. Since (5) implies 
(py)’ S 0 we have even 








102 MATHEMATICS: H. WEYL Proc. N. A. S. 


pe £1,08¢ 81, 


and may therefore start with the set G of all functions g(z) defined and 
continuous in the interval 0 < z < aforwhichO Sg <1. This limitation 
of g has the consequence that 


Osf' s2,0Sf Ss 1/22* (8) 


and enables us to ascertain a universal bound for ¢’, i.e., a number A de- 
pending on \ and a only such that 


0s -y (2) SA 


for 0 S z S a and the image yg = &{g} of everygeG. First we deter- 
mine such a universal bound for 9’(0). Indeed would ¢ start off at s = 0 
with too steep a decline it could not help hitting the ground, y = 0, long 
before z reaches a, considering that the coefficients of its differential equa- 
tion (3) are bounded by (8). Once y’(0) has been penned in that differ- 
ential equation does the rest. Hence the images ¢ of all elements g ¢ G are 
equi-continuous in the strict sense that |z: — 22| < implies |y(z:) — ¢(z)| 
< Ae. 

We topologize the “space”’ of all continuous functions of z by interpreting 
convergence as uniform convergence in the interval0 S z <a. Then # 
is a continuous operator in G which maps G into the subset G,4 consisting 
of all functions g « G for which |z, — %| < ¢ implies lg(a1) — 2(z2)| S Ae. 
This suffices for the existence of a “fixed point” g e G of the operator #, 
(see Birkhoff and Kellogg, l.c., Theorem II). 

The case of the infinite interval is not essentially more difficult. It 
turns out that the solution f of (1) which our method yields has the prop- 
erty that f” = 0 and f’”’ + 2ff” < O converge to zero with z— = at least 
as strongly as a function of the type e~”"(y > 0). Hence 


So°f'(z)dz = f'(o) = B 
converges and we may retrace the step leading from (A),) to (1). 
THEOREM. For given positive k the problem (A,) has a solution w whose 
derivative is monotone increasing from O to k as z travels from O to infinity; 
the second derivative decreases monotonely from 


w"(0) = ak’ (a = p~"”) 


to zero, approaching zero with z — ©& at least as strongly as a function of the 
type e~™. 


1 These PROCEEDINGS, 27, 578-583 (1941). 

? Birkhoff, G. D., and Kellogg, O. P.,. Trans. Am. Math. Soc., 23, 96-115 (1922); 
Schauder, J., Studia Mathematica, 2, 171-180 (1930); Leray, J., and Schauder, J., Ann. 
Sc. Ec. Norm. Sup., 51, 45-78 (1934). 
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UNSTABLE MINIMAL SURFACES WITH ANY RECTIFIABLE 
BOUNDARY 


By Max SHIFFMAN 
COLLEGE OF THE CiTy oF NEw YorK 


Communicated December 17, 1941 


1. Introduction.—Let T be any rectifiable curve in space. We shall 
prove that [ must bound unstable minimal surfaces whenever it bounds 
several minimal surfaces which are relative minima in area. More gener- 
ally, the Morse relations apply in a slightly modified form. In previous 
literature, these results were obtained for special types of rectifiable 
curves.! 

The method is first to obtain unstable minimal surfaces bounded by 
curves approximating I’ and second to perform a passage to the limit. The 
first step may be done in several ways. In the present paper it is carried 
through by extending the procedure used by Courant for polygonal bound- 
aries.» The passage to the limit is accomplished by means of theorem 1 
below, a theorem which was discovered by Morse and Tompkins,* and 
which is here proved by using the isoperimetric inequality for minimal 
surfaces.‘ 

Our method has the advantage of not requiring any explicit representa- 
tion of the Dirichlet integral of a potential surface in terms of its boundary 
values. 

2. Continuity of Area of Minimal Surfaces —2.1. Let r(r, 6) be a po- 
tential surface defined over the unit circle r < 1, with boundary values 
r(1, 6) which are continuous and of bounded variation. Indicate the 
length of the curve z(r, @), for fixed ry and 0 < @ S 27, by L(r); indicate the 
length of r(r, 6) for fixed @and p < r < 1 by L(6; p,1). Thesymbol D[r] 
stands for the Dirichlet integral 


D{t] = a 7 («* + vt dra 


Lemma 1: Li(r) S L(1). 

Lemma 2: Let r"(r, 0), n = 1, 2, ..., be a sequence of potential surfaces 
converging uniformly® to the potential surface r°(r, 0), and suppose that 
L"(i) > L°(1)asn— ~. Then for any « there is a 5 such that any arc con- 
centric to the unit circle whose angle is less than 5 is mapped by each x", n = 
1, 2, ..., ©, into curves of length less than e. 


4D 1\" 
Lemma 3: L(0; p, 1) S (Fu te : log *) except for a set of @’s of total 


measure < 6/2. 
Lemma 1 is proved by use of the Poisson integral. Lemma 2 is a conse- 
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quence of Lemma 1 and the lower semi-continuity of length. Lemma 3 is 
obtained by using well-traversed estimates of D[r].° 
Define a 6-dense radial network in the annular ring p S r < 1 asaset of at 


4 ee : s 
most |= | + 1 radial lines @ = constant, such that two successive radial 


lines have an angular difference of at most 6. Lemma 3 shows that there is a 
5-dense radial network in the annular ring p < r < 1 on each line of which 
the inequality of Lemma 3 holds. 

2.2. We consider a minimal surface r(r, 6) as given in isometric repre- 
sentation, i.e., r is a potential surface and E = G, F = 0 where E = f,?, 
F = 14,°1/r %, G = 1/r? x)”. The area of a minimal surface is equal to its 
Dirichlet integral. 

THEOREM 1.” Let r"(r, 0), n = 1, 2, ..., be a sequence of minimal surfaces 
converging uniformly® to the minimal surface r° (r, 0), and suppose that L”(1) 
—L°(l)jasn— ©. Then 


D[t"] > D[r*] asn > . 


Proof. Let M be an upper bound of the quantities Z"(1), [Z"(1)]*/4z, 
n=1,2,..., ©. By the isoperimetric inequality A < L?/4x for minimal 
surfaces, we have D[r”] < M for all n. 

For any given ¢ choose the 6 of Lemma 2, and then p so near 1 that 


2M (4M 1\'7 4 4M 1 
_ log - + | — + 1) —— -blog - <<. 
T 6 p 6 ua) p 





Consider any x”, 2 = 1,2, ..., ©. For this x”, there is a 5-dense radial 
network in the annular ring p < r < 1 satisfying the inequality of Lemma 3. 
Consider the ith cell of this network. Indicate the Dirichlet integral of r” 
over this cell by D;’; and indicate the lengths of the images under r” of the 
four boundary lines of this cell by a;, c;, 5;, c; — 1, where c;, ¢; — ; correspond 
to the radial lines, a; to the arc of the unit circle and 5; to the are of the 
circle of radius p. Set]; =a; +0;+¢;+¢-,. We havea;, b; < € by 

4M de 
Lemma 2 and the choice of 6, and ¢;,¢; _; < ( - log *) . By the iso- 
p 
perimetric inequality, 


4nD;' <1? = (; +o, +6 -)@+5) + (G+teo¢--1)?s 


1/2 
< E +4 (tog ) | (a; + 5) + — log (1) 
p p 





5 


Denote the Dirichlet integral of r” over the annular ring p < r < 1 by 
D’[{r"]; over the circle r < p by D’[r”]. Summing the inequality (1) for 
at most [47/65] + 1 terms, and making use of 0a; = L"(1) < Mand 0}; = 
L"(p) <= L"(1) < M, we obtain 
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Me . 2M (4M 1\'7 
D'{e"] < = + — (= * log ‘) a 
™ T 6 p 


(* + 1) ie : <(M+1)e. (2) 
6 im) p 


This inequality holds for m = 1, 2, ..., ~. 
Choose NV so large that |D" [r”] — D"[r°]| < eforalln = N. Inequality 
(2) yields 


|\D[x"] — D[x*]| < (2M + 8)e for all n = N. 


The proof of the theorem is completed. 

3. Unstable Minimal Surfaces Bounded by Curves Approximating T.— 
3.1. Let $ be the space of potential surfaces r(r, 6) defined over the unit 
circle, having a finite Dirichlet integral, mapping the circumference mono- 
tonically and continuously onto I, and satisfying the three point condition. 
We define the distance between two surfaces 1, r2 as® 


li — | = Max. = |n(r, 0) — m(r, 0)|. 
rs '/p 


We shall first obtain unstable minimal surfaces in a space @ containing §, 
and then take the limit as § — §. 

Let A = {Aj, Ao, ..., Am} be a collection of m distinct points arranged 
in this order on the oriented curve. Designate the maximum among the 
diameters of the arcs A1A2, AoA;, ..., AmAi1o0f[ by o(A). Thecollection A 
is fixed in this section. 

Let a = {a1, a2, ..., &m_} be a collection of m distinct points arranged 
in this order on the circumference. Define the following three sets of 
potential surfaces: 

(i) B(a) is the set of all surfaces r of $ for which r(a;) = A;,7 = 1,2,..., 
m. 


(ii) R(a) is the set of all convex linear combinations (7, @) of surfaces in 
(a), ie., y is expressible as 


» = ay; where a; > 0, >0a; = 1, and each x; is in B(a). 


(iti) R(a) is the set of all surfaces 3(r, 0) with finite Dirichlet integral ob- 
tainable as limits of sequences of surfaces in R(a), i.e., D[z] < - and 3 = 
lim-y” where each y” is in R(a). The set @(a) may be described as the 
convex closure of the set B(a). 

Concerning §(a), the following statements are easily proved. 

(a) Pla) © Ka) | K(a). 


(b) If 3 is in Q(a), 3(1, 6) is a continuous curve whose length is at most 


the length ZL of [; and 3(1, 6) for a; _; < 6 < a; lies in the convex closure 
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—, 
of the arc A;-1A; of Tf. If rx is any surface in B(a), max. l3(r, 6) — 
t(r, 6)| < ofA). rSl 


(c) The set ®(a) is convex, i.e., if g:, 32 belong to R(a) so does ; = 
(1 — tg + Hi for0 <¢ <1. 

(d) The set R(a) is D-compact, i.e., if 3” is a sequence belonging to R(a) 
for which D[j"] < M, a subsequence can be found which converges to a 
surface 3 in (a). 

(e) Let 3 belong to (a), and vary 3 by setting 


Z'(r, 8) = 3(7, 6) for 0 = 0 + eX(r, 8) 


where X(r, 6) has continuous first derivatives and ¢ is chosen so small that 
led,| < ‘4/2. Transform Z’ by a linear transformation so that the three 
point condition is satisfied, and define 3’ as the potential surface having the 
same boundary values as the resulting surface. Set a’ = image of a under 
these successive transformations. Then 3’ belongs to Q(a’). 

3.2. We are ready to define the space R. § consists of all the surfaces 
in the @(a) for all possible positions of the collection a = {ai, ag,..., 
Gm}. The distance between two surfaces j:, 32 of @ is again defined by 

| — 32 = er laa(r, 6) — 32(r, 6)|. 
By property 3.1 (a), 8 © RK. Property 3.1 (e) shows that it 1s possible to 
perform variations in the space §. 

By using properties 3.1 (c), (d), (e) and methods similar to [3] the fol- 
lowing theorems can be proved: 

THEOREM 2. The problem of minimizing D{3] among all surfaces 3 in 
R(a) has a unique solution, denoted by 3(a). 

THEOREM 3. The surface 3(a) and the quantity d(a) = D{3(a)| depend 
continuously on the collection a = \a,..., &m}- 

As in [3] it may also be proved that d(a) is a differentiable function of the 
Qj, Gg,...5 Ay. 

Define the subspace Q of § as the set of surfaces 3(a) for all possible a. 

Theorem 2 and property 3.1 (c) assert that § may be D-retracted into O, 
i.e., retracted in such a way that the Dirichlet integral never increases. 
Theorem 3 asserts that © is topologically equivalent to the space of all 
possible a’s, and that the Dirichlet integral is continuous in O. 

Thus, so far as Morse theory in R is concerned, we may limit ourselves 
to the space Q; and Morse theory applies to the minimal surfaces in Q. 

In particular, let 11, f2 be any two surfaces in R. The surfaces y;, fe can 
both be D-retracted into surfaces j:, j2in OQ. In O, g and % can be joined 
by a minimizing connected set C; the maximum d of the Dirichlet integral 
on C is the smallest possible among all connected sets in Q joining j:, 32. If 
d > max.{D[rx,], D[re] }, then C contains a minimal surface 3 for which 
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Dj] = d. Because of the isoperimetric inequality for minimal surfaces, 
we have d < L?/4m by 3.1 (8). 

THEOREM 4. Two surfaces t1, 2 in R can be joined by a connected set C of 
surfaces in @ on which the maximum of D{3] is a quantity d < max. {D[xJ, 
D[ro], L?/4r}. If d > max.{D[r], D[re]}, then C contains a minimal 
surface 3 for which D{3] = d. 

4. Passage to the Limit to Yield Unstable Minimal Surfaces Bounded by T. 
—4.1. The considerations in §3 depended on the selection of the collection 
of points A = {A;, As,...,A,} on IT. The space @ defined there will 
now be denoted by @(A). 

Select a sequence A™, A®,..., A™,... of collections of points on I, 
where A” + ” is obtained from A” by adding points of I, and where o(A )) 
—>OQasv— o. Then R(A%) D Q(A®) Dd... DQ(AM)D...DF. 
And using 3.1 (5), lim @(A™) = §. 


4.2. Main Tueorem I. [f the rectifiable closed Jordan curve T bounds 
two minimal surfaces which are proper relative minima, then T must bound an 
unstable minimal surface. 

Proof. Let the two minimal surfaces be 1), %. By Theorem 4, there isa 
connected set C” of surfaces in Q(A™) joining 1, r2 on which the maxi- 
mum of D{3] is a quantity d” satisfying the inequality of Theorem 4. 
Select a sequence of A for which the d converge, and indicate it still by 
A. Consider all the limiting surfaces of the sets C” as vy . They 
form a connected set C in § containing 1, r2 on which the maximum of 
D{tlisd < lim d™ < L*/4r. 


Now, d > max. {D[r,], D[re] } because 1, f2 are proper relative minima. 
Hence d” > max. {D|r,], D[xe] } for all v sufficiently large. By Theorem 
4 there is a minimal surface 3” on C” for which D[;’] = d®. A sub- 
sequence of the 3” converges to a minimal surface 30n C. The lower semi- 
continuity of length and property 3.1 (b) show that the lengths of the 
curves 3’(1, 0) converges to the length Z of [f. Theorem 1 yields 

lim d’ = Dij] < d. 


_o 


Hence d = lim-d” and D{j] = d. The minimal surface required by Main 


Theorem I is 3. 

5. The Morse Relations.—In the same way, by enlarging the space to 
(A) and then passing to the limit 8, one can prove the following: 

(1) any k-cap in § with cap limit a is homologous on $, to a k-cap which 
contains a minimal surface 3 for which D[;] = a;° 

(2) the connectivity numbers of B are RR) = 1,Ri = Re =... =R, = 
... = 0. 

In these two statements, we always consider Vietoris cycles lying on By 
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for some NV. The only difficulty in the proof of these statements, beyond 
the methods already used in $§1-4, is the possibility that surfaces in 
may have boundary values which are constant on some arcs of the circum- 
ference. 

Statement (1) yields a modified Morse theory; the usual requirement is 
that every #-cap contain a minimal surface 3 for which D[3] = a. But by 
assigning type numbers to blocs of minimal surfaces, the usual Morse rela- 
tions hold. 


1 [1] Shiffman, ‘‘The Plateau Problem for Non-Relative Minima,’ Ann. Math., 40, 
834-854 (1940); [2] Morse and Tompkins, ‘‘The Existence of Minimal Surfaces of 
General Critical Type,’”’ Jbid., 40, 443-472 (1940); [3] Courant, “Critical Points and 
Unstable Minimal Surfaces,’’ these PROCEEDINGS, 27, 51-57 (1941); [4] Morse and 
Tompkins, ‘‘Minimal Surfaces Not of Minimum Type by a New Mode of Approxima- 
tion,” Ann. Math., 42, 62-72 (1941); and “The Continuity of the Area of Harmonic 
Surfaces as a Function of the Boundary Representations,”’ Am. Jour. Math., 63, 825-838 
(1941). 

2 Cf. [3]. 

3 Cf. [4]. 

4 See Rado, ‘“‘On the Problem of Plateau,”’ Ergeb. Math., 2, 45-47 (1933). 

5 The uniformity of the convergence may be eliminated since it is a consequence of 
L"(1) > L® (1). 

6 Cf. Shiffman, ‘‘The Plateau Problem for Minimal Surfaces Which Are Relative 
Minima,” Ann. Math., 39, 311-312 (1938). 

7 This theorem is capable of very wide generalizations. 

8 In § this metric is equivalent to the uniform metric. 

9 L, consists of those surfaces I of 3 for which D[Y] < a. 


ON HOMOGENEOUS MEASURE ALGEBRAS 
By DorotHy MAHARAM 
THE INSTITUTE FOR ADVANCED STUDY 
Communicated January 26, 1942 


1. The purpose of this note is to determine the structure of general 
measure algebras. 

2. For any Boolean o-algebra M, let M be the least cardinal number 
which is the power of a o-basis of M. M is homogeneous, if L = M for every 
(principal) ideal L © M different from the null ideal. A Boolean o- 
algebra M = {0 Sa, b,c,... S e} is a measure algebra, if there is defined 
on M a measure function (that is, a countably additive real non-negative 
function) (a) such that (i) 0< ule) < @, (ii) u(a@) = 0 if and only if a = 0. 
We also assume that there is no atomic element in M, i.e., (iii) a > 0 im- 
plies the existence of abe M such that0< b<a 
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An example of a homogeneous measure algebra is the Boolean algebra 
P(y) of all measurable sets (mod. null sets) of an infinite product space 
Q, = Po<a<,J, of intervals [,: 0 S x, S 1, where a and y are ordinal 
numbers and the measure on Q, is defined multiplicatively in terms of the 
ordinary Lebesgue measure on each J,. 

If y < vy’, then P(y) may be considered as a o-subalgebra of P(y’). 
This embedding is obtained by identifying each subset E of Q, with the cyl- 
inder set E X P, <4< +I, in the product space Q,, = 2, X P,<a<yla- 

Two measure algebras are tsomorphic, if there exists a measure-preserving 
o-isomorphism between them. It is easy to see that, for any infinite 
ordinals y and y’, two measure algebras P(y) and P(’) are isomorphic if 
and only if y and y’ correspond to the same cardinal. Moreover, if y is 
finite or corresponds to No, then P(7) is isomorphic to P(1), i.e., to the 
measure algebra of Lebesgue measurable sets (mod. null sets) of the interval 
i: 829 1. 

THEOREM 1. Every homogeneous measure algebra with y(e) = 1 is iso- 
morphic to P(yo), where yo 1s the least ordinal number corresponding to M. 

This theorem will be proved by transfinite induction. It will be suffi- 
cient to prove the following 

Lemma |: Let L bea o-subalgebra of a homogeneous measure algebra M 
with u(e) = 1 such that L < M, and assume that L is isomorphic to P(y), 
where y is an ordinal corresponding to L. Then, for anyae M witha é L, 
there exists a o-subalgebra L' of M such that LCL’, ae L’, and L’ is iso- 
morphic to P(y + 1), in such a way that this isomorphism is an extension of 
the given isomorphism of L and P(y). 

The proof of this lemma will be given in section 4. 

3. Let us assume the conditions of Lemma 1. We can consider Q, 
as a representation space of L = P(y). Then, by a theorem of Radon- 
Nikodym, there exists, for any a e M, a measurable function ¢,(¢) defined 
on Q, such that 0 S ¢,(§) S 1, and Si,a(é)dé = p(a A x) for any xe L, 
where E, is a measurable set of 2, which corresponds to x e L. An element 
ae M is called independent of L, if we have ¢,() = constant [namely, = 
u(a)] almost everywhere (a. e.) on 2,. This is equivalent to saying that 
we have u(a A x) = w(a)-u(x) for all x « L. 

Lemma 2: For any measurable function x(&) defined on ©, with 0 S x(é) 
< ¢,(£) a. e. on Q,, there exists abe M such that0 S 6 Saand ¢,(E) = x(§) 
a. €. on Q,. 

Proof of Lemma 2.—Our lemma is clear if x(§) = 0 a.e. on Q,. Hence 
we may assume that meas[é:x(é) > 0] > 0. Because of the principle of 
exhaustion, it is sufficient to show that there exists a )* « M such that 0 < 
b* < aand g+(é) S x(é) a.e. onQ,. 

Let A be the principal ideal generated by a and let L(a) be the o-subalge- 
bra of M generated by Landa. Since L(a) = L< A = M by assumption, 
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there exists a b; eM such that b, e A and },¢ L(a). This means that meas[é: 
0 < ,(&) < ¢,(&] > 0. Again, by the principle of exhaustion, we can 
further find a b, e M such that 0< b, < a and0< %,(&) < ¢,(¢) a. e. on the 
set [E:9,(£) > O]. Let us denote by c, and @ the elements of L which 
correspond to the sets [£:0 < g,(£) S 27'y,(£)] and [€:2~*y,(£) < o,(é) < 
¢a(€)]. If we now put bs = (4, A be) V (ce A (a — be)), then0<d<a 
and we have 0 < ¢;,() S 2~*¢,(£) a. e. on the set [E:¢,(£) > 0]. By iter- 
ating this argument, we can find, for each n, a b, +2 ¢ M such that 0 < d, 4.2 
<aand0< %,,,() S 2~"o,(£) a. e. on the set [E:9,(£) > O]. Taken 
so large that meas E* > 0, where E* = [£:2~",(&) S x(&)], and denote by 
c* the element of L which corresponds to E*. If we put b* = c* A dy 4., 
then we have 0 < b* < a, and g(t) = by + 26) Ss 2~"o,() S x(é) a.e. 
on E*, and g«(€) = 0a.e.on Q, — E*. Hence gys(~) < x(é) a. e. on 
Q, and this proves Lemma 2. 


4. Proof of Lemmal. Let A, be the set of all finite sequences 6 = {a1, 
..+> €,}, Where e; = Oor 1,7 = 1,...,m; and let A = U,P_,A,. A count- 
able set {as} (6 « A) is a dyadic decomposition of a if a V a =a, 
do A a, = 0, and ag 9 V 45,1 = 45,45,9 A a3, = Oforallée A. By Lemma 
2, there exists for any a « M, a dyadic decomposition {a;}(6 « A) of a 
such that ¢,,(£) = min.(¢,(£), 6/2 + ... + «,/2" + 1/2”) — min. (¢,(6), 
4/2 +... + ,/2")a.e.onQ, for alld = {a,...,¢,}€A. In the same way, 
there exists a dyadic decomposition {a’,}(6 « A) of a’ = e — a such that 
$o',(€) = max. (¢a(€), «,/2 ee €,/2” + 1/2”) — max. (go(€), «,/2 + 
... + ,/2”) a.e.onQ; for all 6 = {a,..., &}eA. Let us put b; = a; V 
a’; for alld ¢ A. Then {b;}(6€ A) isa dyadic decomposition of the unit 
element e, and each ), is independent of L, since we have clearly ¢,,(¢) = 
$a,(€) + ¢a’,(€) = 1/2” a. e. on 2 for all d¢ Ay, m = 1, 2,.... 

Now let L’ be the o-subalgebra of M generated by LZ and (b,}(5e A). 
L’ is obtained by completing the finite algebra L* consisting of all elements 
of M of the form: b = V5ea,(bs A ), where ce L for all 6 A,,m = 
1, 2.... It is easy to see that L’ is isomorphic to P(y + 1) by an iso- 
morphism which is an extension of the given isomorphism of L and P(7). 

Finally, in order to prove that a e L’, consider the set [:9,(€) 2 4/2 + 
... + €,/2" + 1/2"] for each 6 = {a,..., &,}€A, and let cs be the corre- 
sponding element of L. If we putd™ = v4, 4, (0s A ¢), then we have bY < 
b <... 5b <<... Saand ula — b™) S 1/2" forn = 1,2, .... Hence 
we havea = v2. ,b™ and consequently ae L’. This proves Lemma 1 
and so Theorem 1. 

5. A measure algebra M is a direct sum of (a finite or countably infinite 
number of) measure algebras M,, if each M, is (isomorphic to) a principal 
ideal of M and if every element a « M is uniquely expressed in the form: 
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A= Van w= On Im A Ay = O(m ~ n), wherea,¢M,,n = 1,2,.... (Thislast 
condition is equivalent to saying that the unit elements e, of M,,, which are 
elements of M, satisfy ey, A €, = 0(m An) ande = V pa i€p-) 

THEOREM 2. Every measure algebra is a direct sum of homogeneous 
measure algebras M,,(n = 1, 2, ..., finite or countably infinite). 

The proof of Theorem 2 is easy and will be omitted. 

Theorems 1 and 2 indicate the structure of a general measure algebra. 

6. The ergodicity of a measure preserving o-automorphism T of a 
measure algebra M (onto itself), or that of a group G = {T} of such o- 
automorphisms, can be defined as usual. 

Lemma 3. In order that the group of all measure preserving o-automor- 
phisms of a measure algebra M be ergodic on M, it is necessary and sufficient 
that M be homogeneous. 

The proof of this lemma is easy and is omitted. 

THEOREM 3. Let M be a measure algebra, and let G be the group of all 
measure preserving o-automorphisms of M. Then M is a direct sum of a 
countable number of invariant principal ideals M,, on each of which G its 
ergodic. This decomposition is the same as in Theorem 2. 

Let G be an arbitrary group of measure preserving c-automorphisms of 
M. Then the set L¢ of all a « M such that T(a) = a for all T € G, isa o- 
subalgebra of M. Conversely, for any o-subalgebra L © M, consider the 
set G, of all measure preserving o-automorphisms T of M such that T(a) = 
aforallaeZ. Gz, isclearly a group. It is also clear that L © Lg, and 
G ©G_,,. Exactly as in the theory of Galois, we may ask the question: 
Under what conditions do the equalities L = Lg, and G = G,,, hold? 
These and related problems will be discussed elsewhere. 

I am greatly indebted’ to Shizuo Kakutani, who revised and improved 
my original manuscript. 
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ON THE DEPOLARIZATION OF NEUTRON BEAMS BY 
MAGNETIC FIELDS 


By O. HALPERN AND T. HOLSTEIN 


RADIATION LABORATORY, MASSACHUSETTS INSTITUTE OF TECHNOLOGY AND WESTING- 
HOUSE RESEARCH LABORATORY 


Communicated January 10, 1941 


In two previous notes! referred to henceforth as I and II the authors have 
discussed the transmission of neutrons through magnetized and unmag- 
netized ferromagnets and have in detail reported on the physical conse- 
quences derivable from observations on polarized neutron beams. For all 
the applications reference should be made to II. 

The mathematical method used was semi-classical; the positional co- 
ordinates r of the neutron were treated fully classically because the wave- 
length is very small compared to the linear dimensions of the domains in 
which the magnetic field changes; the behavior of the spin was described 
with the aid of the self-explanatory equation 


= = gis x BO)! (1) 

t 

which could be integrated for all physically interesting cases as far as the 
dependence of B on r was concerned; the result of this integration, when 
averaged over all neutron paths, gave directly the change in the average 
spin. 

There exist problems of considerable physical interest in which it is 
necessary to know not only the average spin but the average number of 
neutrons in each individual spin state. For this type of problem the 
method described and applied below seems to be appropriate. We shall 
limit the discussion to two cases which are (cf. II) of direct physical signifi- 
cance. The method is general and applies to all kinds of particles and all 
values of the spin. 

The change in the number of neutrons occupying each spin state is deter- 
mined as follows: Every individual particle undergoes changes in its spin 
state under the influence of the magnetic fields. The final result is ob- 
tained by averaging over all paths of the various particles which in general 
will traverse regions of different magnetic field. It will be seen below that 
the effect of the transmission is to equalize the average number of neutrons 
in each spin state; since a polarized beam is characterized by a preferential 
occupation of the different spin states, the effect of a ferromagnetic medium 
is depolarization of such a beam. 

(1) The magnetic field B; of the 7th domain which is traversed by a 
neutron during the time interval 7; is random in direction but constant in 
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magnitude and direction during 7;. The times 7; are in addition limited 
by the condition gB;jr; << 1. This case constitutes an excellent approxi- 
mation to the behavior of an unmagnetized polycrystalline specimen, 
whose crystal grains have linear dimensions < < 10-* cm. 

The starting point is the wave equation 


dc, 


th ad = LA amen: (2) 


Here c,, is the amplitude of the state for which s, = nh and the H,,,, are the 
matrix elements of the Hamiltonian gs-B viz.: 


Hyun = nhgB, (3a) 





He gai ™ 5 gh(B, 7 iB Viz nltn+i)-. (30) 


where 


_ {se|max. 


l 
h 


As before, the positional codrdinates are treated classically : 
x = X% + vt, y = Yo, 2 = %. 


The existence of randomly oriented fields is taken into account in the fol- 
lowing manner: The total time 7 is divided into intervals 7; each of 
which is characterized by a Hamiltonian matrix H®, the average of which 
over all neutron paths has the following properties: 


H®, = 0, HEHE = thu. (4) 


nm7i 


Hi, 
The condition gB;r; << 1 becomes —— qe hs 
i 


One could now find exact solutions of (2), and average by (4); this pro- 
cedure is somewhat cumbersome and can be avoided by the following 
method. One introduces a time interval 3} intermediate in magnitude be- 
tween the 7; and TJ such that: (@) many domains are traversed in #; 
(b) the relative change suffered by c, during 3? is small. Using the second 
property, one can immediately write down the values of the c, at a time 
t + 3 by a method of successive approximations: 


c(t + 8) = A(t + 8) + ME + 8) + Pe + ) (5) 
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where 
Cult + 3) = c(t) 5a) 
Ot + 8) = 5 LD Hem (tr (56) 
t+ aA=-5 a 2 2d HH Cy(t)ti75 — 


h? i j<im 
> . o 
jp PLL HamHmace(l)ri. (5e) 


> here indicates summation over all intervals 7; which are enclosed in #. 


$ 
The first term of (5c) represents the effect in domains (7) of the changes in 
C, produced in previously traversed domains; the second term gives the 
second order effect of the ith domain, itself. For the probability w,(¢ + 8) 
= ¢,c% one obtains, to the second order: 


Wa(t + 8) = w,(t) + [co"(t)c (t + 8) + comp. conj.] — 


he < z HamE pm’ (t)Cm? (t) 747; arg 


hh? i m, m’ 


1 ‘ : 
i b DY AYA Yc, (t)ex(t) + compd. con. ~- 
i m,? 


bl z. > HYHYc,(t)ck(t)7;7; + compl. conj. |. (6) 
2H? 4, j<i m, 4 
For the average over all neutron paths, i.e, for the average over all directions 
of the magnetic field, one obtains by (4) 

] s ’ 
y(t + 8) = wy) + =D HGH Goenlt)cne(rt — 


1 


oye b LH aml mpcy(t)en(t)7 t; + compl. conj. |- 


wy(t) + 5 , DSP Wt — 5 > ASP wae (7) 


In obtaining the last term, the Hermitian property Hy», = Hii, was 
employed. 

Since # is so chosen that only small changes occur in c, and, consequently 
in w, one can replace the differences by differentials, obtaining: 





ry 
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) 2 —_—. pees 

> a z i. D | Hnm|?(m — Wy) = = Lol Anm|?(m — Wy). (8) 

t h? +; ‘m hi? ‘ns 
This equation is an example of a type commonly encountered in statistics. 
The change in the probability of occupation of a state, m, is equal to the 
probability of transitions from all other states, m, minus the probability of 
reverse transitions—from tom. A stationary condition is reached when 
all the states are equally occupied. The manner in which equilibrium is 
attained is somewhat complicated for the general case, but can be illus- 
trated very simply for the particular case of spin '/2. Equation (3) to- 


gether with the equality of B?, B} B? leads to the equations: 


ow, 1 


1 Ow) 
-* = f 


ot 





(9) 


*7Bw_s), — Wy) = — 


The connection with the classical method can be established by calculat- 
ing the expectation value for the polarization. In the simple case of spin 
1/2 this polarization s, is s/2(w:;, — w_1,,); the equation for its variation is 
easily obtained from (9): 


— = — - g*rB;s, (10) 
which has the solution 
yf pone 
S= Sqo EXP (- 3 e'rBit) (10a) 


to be compared with II. 
The equation for the spin can readily be obtained in the general case of 
spin/. Through use of (3) equation (7) becomes 


ow, 1 


eae 72°(Bz + By)[(l — n)(L + 0 + 1) (We 41 — Wn) + 


The equation for the z component of the spin is: 


Os, re) 1 
=A +. Lnw, = qerB: + Bs) Vi [nw, 4 (2 — n)\(-+n +1) + 


NW, —~,(1+n)(l—n+1) — nw, {(l —n)(+n+1) + 


7+n)¢i—n+1)}]. (12) 
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Now 
n’=1+1 
Yinw, +10) -—1)(l+n24+1) = > . (n' —1)—n'+1)wy = 
n=l 
yy An —-1)\¢-—n+1)¢+n)u, 


since the terms not common to both summations actually have the value 
zero, as can be learned from inspection. Similarly 
n=l ntl 


> nv,-il + nl—-n+1)= DY w+) —nlt+n4 l)y,. 


n=—1 n=—l 
Substituting in (12) one obtains 


Os, 1 — — 
= = oe + eo = — 


i 7g°B?s, (13) 


wie 


which has the solution 
1 toes 
S,(t) = Sy Exp (- 3 r¢'B*t) (13a) 


to be compared with (10a). 

This formula constitutes the quantum-mechanical verification of the 
important fact that the depolarization depends, not on the magnetic moment 
or the spin separately, but only on their ratio g. 

(2). The total field B of each domain is the sum of two fields By + B;. 
The first, Bo, is constant throughout the time 7. The second B; is con- 
stant within 7; The direction of B; is perpendicular to that of By but 
otherwise random while the absolute value of B; is small compared to that 
of Bo. It will furthermore be assumed that the initial direction of the spin 
is parallel or anti-parallel to By. These assumptions (cf. II) approximate 
the conditions prevailing during the passage of neutrons through a quasi 
saturated polycrystal. 

These conditions take on the following mathematical form: 





B,; = 0, B,,By; = 0 unless i = j, x = y 
=e — Pa re ae Fy — 
5. bic B,; ia 0; B; fy By; or 5 Bi. 

Due to the fact that gBr; is not small, an approach analogous to that of 
the preceding section is not applicable. In order to solve (2) explicitly, one 
sets up a time-dependent perturbation treatment. Introducing the 
quantities a, = c, exp (—iE,t/h) where E, = H,,, one obtains the well- 
known equations 








Ci 
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. da, . 
th 7 DV Aumomexp[—-t(E, —  ,)t/h). (14) 
If the axis of quantization is chosen parallel to the vector Bo 
H,, = E, = nghBo (15a) 
He, 21 = Yehg(By; = iBy) V (1 = n)(l=nm+1) (158) 


from which it is apparent that the usual mechanism of the time dependent 
perturbation theory is an expansion in the parameter |B;|/|Bo|. 
The solution of (14) is similarly to (5): 


dy(t + 3) = a(t) + aP(t + 8) + oP (t + 9) (16) 
a(t + 8) = OY AMa,, (1) exp (—i(E, — E,)0;/h) 


$ 





exp(—12(E,, rah En)t;/h) molt 














eg ae 
(2) i () 77 @) exp(—i(E, — E,)0/h) 
ay (t + #) ys 2 4 Laml mp p(t) E; eae x. 
ea ae E,)t;/h) uae. nah exp(—7(E, ig E,,)1;/h) “3¥ 1] of 
E, isin (om ff. na Pm | 
2D, LD HamHrag(0) exp[—t/M (Em — En) 9: — (Ey — Bn)03}1. 
EGF <SE Mh, 
exp(—i/B(By ~ En)ri) ~ 1 exp(—i/MEm ~ Ep) 4) 1 gy 
E, — Em En — E» 
j = i-l 
with 0; = bw Tj. 
j=l 


One now sets up the probability, w,(¢), as before, and averages over the 
direction of the magnetic fields B;, obtaining 


a \Hnm|? ae. 
w,(t + 3) = w,(t) + 4 2 y G, — E) (Wm — Wy) sin’ (E, — Em)r;/2h. 
(17a) 


Replacement of finite differences by differentials yields, finally: 
dw, Hy 2 4 

ae OP per eg ee _ nz —s E : , 
a~% (E, — Em)? 7 (Wm — Wy)sin?(E, — En)r;/2h. (170) 





Equation (17b) can be integrated without particular difficulty for any 
case of actual interest. One can on the other hand learn from it without 
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calculation how the two qualitatively different results of II were arrived at. 
We found that for gBr < < 1 the depolarization coefficient is proportional 
to Bj independent of By while in the other case it turned out to be propor- 
tional to Bj/Bj. These results are here immediately obtained from 
(17) since it is legitimate in the first case to expand the sin? thereby remov- 
ing the energy denominator (E, — E,,)? which is proportional to Bj. In 
the second case sin? may be replaced by 1/2 leaving thereby the energy 
denominator ~ B;”. 

These mathematical features permit an informative physical interpreta- 
tion. In each time interval 7; matrix elements ~gB,;i cause transitions 
between energy levels which are separated by ~gfiBy. The uncertainty 
principle AEr ~ i when applied to this case states that transitions occur 
with appreciable frequency only between such states which permit con- 
servation to within i/r. If therefore gB,i > h/r then transitions are cut 
down; which leads mathematically to a polarization coefficient propor- 
tional to B}/Bj. If, however, gBosh << ii/r, then no restrictions on the 
frequency of transitions are imposed by the energy principle since the 
quantum theoretical uncertainty exceeds the energy differences between 
the various states. The frequency of transitions is then solely determined 
by the values of the non-diagonal matrix elements which are proportional to 
B; and independent of Bp. 


1 Halpern, O., and Holstein, T., (I) Phys. Rev., 55, 601 (1939); (II) Jbid., 59, 960 
(1941). 








